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A  bs  tract 

"A  method  is  developed  for  obtaining  uniformly  valid  so¬ 
lutions  about  slender  bodies  (2-D  or  j-U)  with  round  lead¬ 
ing/trailing  edges  in  compressible  subsonic  flow.  The  meth¬ 
od  is  based  on  the  principles  of  singular  perturbation  theo¬ 
ry  and  applied  to  nonlifting  thickness  problems.  The  method 
corrects  small  perturbation  theory  (outer  solution)  with  a 
stagnation  region  theory  (inner  solution)  yielding  uniformly 
valid  solutions.  Analytic  results  of  the  first  and  second 
order  method  are  presented  and  compared  to  other  methods, 
including  exact  numerical  results,  for  the  elliptic  cylinder. 
The  first  order  method  is  shown  to  be  preferred  to  higher  or¬ 
der  efforts  and  is  directly  applicable  to  current  numerical 
techniques  solving  the  Prandtl-Glauert  equation.  Application 
of  the  first  order  method  numerically,  using  a  current  panel 
technique,  is  presented  for  the  NACA  0012  airfoil  and  a  glide 
bomb  configuration.  All  comparisons  of  the  first  order  com¬ 
posite  results  with  other  methods  sho w  excellent  agreement 
and  improved  results  subject  to  the  limits  of  validity  of  the 
theory.  The  ability  to  correctly  model  the  flow  behavior  in 
stagnation  regions  using  the  composite  solutions  leads  to  a 
greater  confidence  in  aerodynamic  predictions,  especially  in- 
viscid  drag  estimates. 
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The  purpose  of  this  dissertation  is  to  develop  a  method, 
based  on  the  principles  of  singular  perturbation  theory,  for 
obtaining  uniformly  valid  solutions  for  flow  about  slender 
bodies  with  round  leading  or  trailing  edges  in  subsonic  po¬ 
tential  flow.  Specifically  the  method  improves  solutions 
obtained  from  the  linearized  small  perturbation  potential 
equation  (Prandtl-Glauert  equation)  in  stagnation  regions. 

The  ability  to  correctly  model  the  behavior  in  stagnation 
regions  leads  to  a  greater  confidence  in  aerodynamic  predic¬ 
tions,  especially  inviscid  drag  estimates.  The  first  order 
method  is  shown  to  be  preferable  to  a  second  order  effort, 
and  is  directly  applicable  to  current  numerical  techniques 
solving  the  Prandtl-Glauert  equation.  Results  of  the  first 
order  method  are  compared  to  numerical  solutions  of  the  exact 
potential  equation,  and  are  shown  to  be  in  excellent  agree¬ 
ment  subject  to  the  limits  of  validity  of  the  theory.  The 
present  development  is  for  nonlifting  bodies  only. 


Background 

In  recent  years  there  has  been  considerable  theoretical 
and  engineering  development  of  numerical  techniques  for  the 
prediction  of  potential  flow  about  arbitrary  configurations. 
For  compressible  flows  the  most  useful  methods  to  date  for 
arbitrary  configurations  (2-D  or  3-D)  have  been  the  surface 
singularity  methods,  or  panel  methods  (Ref  1-6).  The  solu¬ 
tions  obtained  by  these  methods  for  steady  compressible  flow 
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are  solutions  to  the  well  known  1‘rmni tl-Glauert  equation. 

'Phi a  oquntion  i : i  dorivod  from  l.ho  ux.irl.  potential  equation 
based  on  the  assumption  of  small  perturbations  or  distur¬ 
bances,  and  is  valid  for  subsonic  or  supersonic  flow  (Ref  7). 
Solutions  thus  obtained  offer  reasonable  predictions  of  in- 
viscid  lift  and  pitching  moment,  and  surface  pressure  dis¬ 
tributions  except  in  stagnation  regions  (Ref  8,9).  The  fail¬ 
ure  of  these  methods  in  stagnation  regions  often  results  in 
the  prediction  of  unrealistically  high  surface  velocities, 
and  thus  poor  surface  pressure  estimates  in  these  regions. 

The  overall  result  is  poor  inviscid  drag  estimates,  espe¬ 
cially  for  configurations  with  round  or  blunt  edges,  where 
the  local  slopes  are  greatest  (drag  being  the  integral  of 
pressure  times  slope  in  the  free  stream  direction).  This 
failure  is  directly  related  to  the  fact  that  the  assumption 
of  small  perturbations  in  the  stagnation  region  is  very  poor. 

History  of  Previous  Efforts 

Thin  airfoil  and  slender  body  theory  provided  aerodynam- 
icists  the  first  opportunities  to  make  accurate  estimates  of 
the  potential  flowfield  about  arbitrary  conf igurations ,  espe¬ 
cially  for  those  immersed  in  a  compressible  flow.  It  was 
recognized  by  many  investigators  however,  that  these  results 
were  not  accurate  in  stagnation  regions.  A  method  to  correct 
thin  airfoil  theory  in  stagnation  regions  was  developed  by 
Lighthill  (Ref  10),  while  an  alternate  approach  obtaining  the 
same  results  was  developed  by  Van  Dyke  (Ref  11).  These  meth¬ 
ods  are  typified  by  Riegel’s  rule  (Ref  ll').  With  the  advent 
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of  electronic  computers,  and  the  more  exact  modeling  obtain¬ 
able  with  surface  singularity  methods,  the  above  mentioned 
methods  are  rarely  used.  Current  emphasis  seems  to  be  di¬ 
rected  at  correcting  the  resultant  inviscid  drag  estimates 
based  on  the  theoretical  results  of  Jones  (Ref  12,13)  as 
presented  by  Sotomayer  and  Weeks  (Ref  14).  This  approach 
involves  estimating  the  amount  of  leading  edge  suction,  and 
thus  leading  edge  thrust  or  drag,  necessary  to  obtain  accu¬ 
rate  results.  A  much  different  and  fortuitous  approach  de¬ 
veloped  by  Stancil  (Ref  8)  for  improving  supersonic  wave 
drag  estimates  could  also  be  used  for  subsonic  applications. 
This  method,  deemed  modified  linear  theory,  solves  the  Prandtl- 
Glauert  equation  based  on  the  local  Mach  number  instead  of 
the  free  stream  Mach  number.  Each  of  these  methods  will  now 
be  discussed  in  more  detail,  with  acknowledgements  to  each 
for  providing  insights  into  the  present  method. 

As  mentioned,  the  methods  developed  by  Lighthill  and  Van 
Dyke  are  for  correcting  the  results  obtained  from  thin  air¬ 
foil  theory  in  stagnation  regions.  Since  Lighthill* s  method 
is  valid  only  for  round  edges,  and  his  results  are  duplicated 
by  Van  Dyke,  only  the  method  as  presented  by  Van  Dyke  will 
be  discussed.  Van  Dyke's  method  is  based  on  the  principles 
of  singular  perturbation  theory  and  the  theory  of  matched 
asymptotic  expansions.  Two  regions  in  the  flowfiald  about 
a  slender  body  are  clearly  identified,  the  dominant  region 
being  one  of  small  disturbances  from  the  free  stream  condi¬ 
tions  and  the  other  being  the  stagnation  region.  The  results 
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of  thin  air-loj  l  theory  aro  salvaged  (Kef  11)  by  applying  a 
multiplicative  corredtion  which  is  dependent  on  the  shape  of 
the  surface  in  the  stagnation  region.  For  round  edges  in 
incompressible  2-D  flows,  the  flowfield  in  the  stagnation 
region  is  approximated  by  the  exact  incompressible  flow  about 
a  parabola,  chosen  to  match  as  closely  as  possible  the  given 
geometry  in  the  stagnation  region.  This  forms  the  basis  for 
the  correction.  For  round  edges  in  compressible  flows,  the 
exact  solution  for  the  flow  about  a  parabola  is  not  known 
and  must  be  approximated  using  a  Janzen-Rayleigh  expansion 
or  experimental  data.  The  thin  airfoil  results  themself  are 
of  course  extended  to  compressible  flow  using  the  Prandtl- 
Glauert  compressibility  correction.  This  method  for  sharp 
edges  follows  the  same  approach  using  the  flow  about  a  wedge 
to  form  the  basis  for  the  correction.  The  region  of  non¬ 
uniformity  for  sharp  edges  in  compressible  flow  is  so  small 
that  any  correction  is  negligible  and  therefore  not  used. 

Thus  the  method  just  described  is  applicable  only  to  the  re¬ 
sults  obtained  from  thin  airfoil  theory  (2-D),  and  exact 
only  for  incompressible  flow.  Extensions  to  the  more  exact 
modeling  available  using  surface  singularity  methods,  and  to 
3-D  configurations,  has  not  been  adequately  developed. 

The  method  developed  by  Sotomayer  and  Weeks  to  improve 
inviscid  drag  estimates  obtained  using  panel  methods  is  rep¬ 
resentative  of  methods  based  on  the  theoretical  results  of 
Jones.  Jones  showed  that  when  applying  thin  airfoil  theory 
to  an  airfoil  with  a  rounded  nose  in  incompressible  flow. 


4 


the  resultant  pressure  distribution  in  the  stagnation  region 
is  in  error,  in  evaluating  the  drag  using  those  results  one 
mu3t  add  a  leading  edge  drag  term  to  obtain  the  correct  in- 
viscid  2-1)  result,  that  boing  zero.  (Note  that  for  the  lift¬ 
ing  problem,  or  flat  plate  at  angle  of  attack,  one  must  add 
a  leading  edge  thrust  term. )  The  error  is  due  to  a  leading 
edge  singularity  which  is  an  essential  part  of  the  thin  air¬ 
foil  theory  solution.  The  method  of  Sotomayer  and  weeks 
attempts  to  extend  the  correction  (leading  edge  thrust/drag) 
developed  by  Jones  for  2-D  incompressible  flow  to  J-D  wings, 
correcting  for  compressibility  and  sweep.  For  2-D  flows  the 
amount  of  leading  edge  thrust/drag  required  is  easily  deter¬ 
mined  since  one  knows  beforehand  the  desired  result  (zero), 
but  for  3-D  flows  (induced  drag)  this  is  not  the  case  and 
one  requires  some  other  criteria  to  base  the  correction  on. 
Normally  other  empirical  or  experimental  data  for  similar 
configurations  is  used  as  this  basis  for  determining  the 
level  of  the  correction  required.  Thus  while  this  technique 
can  produce  favorable  results  in  comparison  with  test  data, 
the  difficulty  arises  in  predicting  the  amount  of  correction 
necessary.  Although  there  are  guidelines,  the  amount  of 
leading  edge  thrust  necessary  can  vary  from  zero  (po  correc¬ 
tion)  to  100%  of  the  theoretical  value,  and  is  by  no  means 
analytically  determined. 

Stancil's  modified  linear  theory  was  developed  to  im¬ 
prove  supersonic  wave  drag  calculations,  but  could  also  be 
applied  to  subsonic  flows.  The  foundation  of  the  method  is 


in  the  brand  Li-(Jiuuor t  equation  bailed  on  iocai  Mach  number 
instead  of  the  free  stream  Mach  number t 

(  (  -  M  ’ )  (px<  *-  cpY v  +  (f)„  0 

On  examining  this  equation  it  is  evident  that  as  the  local 
Mach  number  (or  velocity)  approaches  zero  in  the  stagnation 
region,  Laplace's  equation  results.  On  the  other  hand  for 
flow  ovor  a  slender  body,  whore  the  local  Mach  number  is 
approximately  the  free  stream  value,  the  Prandtl-Glauert 
equation  results.  (This  again  points  out  the  diverse  na¬ 
ture  of  the  various  flowfield  regions. )  The  method  is  ap¬ 
plied  using  current  panel  techniques  with  surface  singular¬ 
ity  distributions.  The  exact  surface  velocity  boundary  con¬ 
dition  (as  opposed  to  thin  airfoil  type  approximations)  is 
applied  to  help  eliminate  unrealistic  negative  peak  pressure 
near  leading  edges.  The  method  requires  an  iterative  scheme 
and  yields  attractive  results.  There  are  discrepancies  how¬ 
ever  in  the  solution  method  especially  using  panel  or  sur¬ 
face  singularity  distributions  which  are  a  function  of  ^3 
(l/^3  to  be  more  exact).  First,  when  determining  the  effects 
of  other  panels  on  a  particular  panel  (control  point),  the 
local  fi  for  that  particular  control  point  is  used,  making 
each  panel  solution  at  best  a  local  solution.  Putting  these 
solutions  t-.ogether  globally  is  questionable  and  has  not  boon 
rationally  justified  by  the  author.  Second,  as  the  Mach 
number  approaches  one  (which  must  occur  somewhere  on  the 
surface  for  supersonic  flows),  l/fi  approaches  infinity.  To 
avoid  the  obvious  difficulties  this  wouldcause,  the  author 
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in Lroducoti  a  correlated  Local  Mach  or  p  { baaed  on  calcula¬ 
tions  lor  L!-U  ramp:;  and  cones)  Lor  l.heso  calculations.  Thus 
while  the  method  yield:;  improved  results,  it  does  not  rest 
on  the  nound  theoretical  or  rational  foundation  desired  be¬ 
fore  general  application. 

Objective  and  Approach 

The  characteristic  feature  of  all  the  aforementioned 
methods  is  that  two  distinctly  different  flow  regions  in  the 
flowfield  about  a  slender  body  are  clearly  evident.  One  is 
a  region  where  the  perturbations  from  the  free  stream  condi¬ 
tions  are  small  as  assumed  in  classic  thin  airfoil  or  slender 
body  theory.  The  other  region  is  the  stagnation  region  where 
these  small  perturbation  assumptions  are  invalid  and  produce 
erroneous  results  there.  The  objective  of  this  effort  is  to 
develop  a  method  for  obtaining  uniformly  solutions,  appli¬ 
cable  to  arbitrary  configurations  with  round  edges  in  sub¬ 
sonic  compressible  flow.  This  of  course  means  that  the  solu¬ 
tion  must  combine  the  correct  behavior  in  the  stagnation  re¬ 
gion  with  that  for  the  remainder  of  the  body  while  preserving 
the  distinctive  nature  of  each. 

Chapter  II  presents  the  theoretical  development  of  this 
new  method.  The  approach  is  guided  by  singular  perturbation 
theory  in  defining  these  two  regions  and  identifying  the 
appropriate  governing  equations  for  each.  First  and  second 
order  problems  are  formulated  including  appropriate  boundary 
conditions.  The  formation  of  a  uniformly  valid  composite 
solution  synergistically  combining  the  effects  of  each  region 
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is  then  proposed  based  on  hypotheses,  developed  in  that  chap¬ 
ter.  Being  essentially  a  per  l.urbri  1. 1  on  procedure,  the  result¬ 
ing  method  is  an  approximate  method  which  relies  on  physical 
reasoning  for  understanding  many  of  the  steps  taken  in  the 
development.  This  is  necessarily  a  consequence  of  the  prob¬ 
lem  at  hand  as  the  exact  potential  equation  is  highly  non¬ 
linear  and  adopting  some  sort  of  perturbation  procedure  ap¬ 
pears  to  be  the  most  attractive  approach  to  obtaining  rea¬ 
sonable  results.  This  is  not  to  say  mathematical  rigor  is 
unnecessary,  for  it  will  be  provided  where  possible  in  this 
development,  however*  “...the  proper  questions  to  ask  and 
approximations  to  make  often  come  from  physical  arguments. " , 
(Ref  15). 

The  application  of  the  method  is  to  elliptic  cylinders 
(2-D)  of  varying  thicknesses  in  subsonic  compressible  flow. 

The  ellipse  was  chosen  to  test  the  approach  with  analytic 
as  opposed  to  numerical  solutions.  Chapter  III  presents  the 
first  order  analytic  solution  (Appendix  A),  application  and 
results.  The  first  order  analytic  results  are  compared  with 
the  results  from  thin  airfoil  theory  (Appendix  B)  and  a  nu¬ 
merical  solution  to  the  exact  full  potential  equation  obtain¬ 
ed  using  a  Jameson  code  (Hef  lb, 17).  The  first  order  analytic 
solution  (the  first  such  solution  with  the  stated  boundary 
conditions)  is  essentially  the  solution  to  the  Prandtl-Glauert 
equation  with  the  mass  flux  surface  boundary  condition  (Ref 
18).  (A  discussion  of  the  mass  flux  surface  boundary  condi¬ 
tion  and  its  necessity  is  presented  in  both  Chapters  II  and 
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Hi.)  The  first  order  results  arc  shown  to  compare  very  fa- 
vornblv  wi  th  the  exact;  numerical  mail  Is  for  thin  bodies  in 
subsonic  compressible  f’Low. 

In  Chapter  IV  the  first  order  method  is  developed  and 
applied  numerically  using  a  current  paneling  program  as  de¬ 
veloped  by  Woodward  (Ref  2),  USSAERO,  version  B.  The  purpose 
of  this  application  is  to  show  that  the  method  can  be  used 
with  any  numerical  technique  that  essentially  solves  the 
Prandtl-Glauert  equation  with  the  mass  flux  surface  boundary 
condition,  and  is  not  restricted  to  a  specific  program.  The 
first  application  is  to  a  NACA  0012  airfoil  (2-D)  and  the 
results  are  compared  to  the  standard  paneling  results,  thin 
airfoil  theory  (Appendix  C),  Jameson  results,  and  experimen- 
ta 1  data.  The  second  application  is  to  a  3-D  glide  bomb  con¬ 
figuration  with  the  results  being  compared  to  the  standard 
paneling  results.  Calculated  inviscid  drag  is  consistently 
lower  than  that  calculated  by  standard  methods.  The  new  re¬ 
sults  are  thus  in  better  agreement  with  the  theoretically 
exact  value  of  zero  to  be  expected  in  subsonic  flow.  This 
leads  to  a  higher  level  of  confidence  in  predicting  inviscid 
drag  levels  for  complex  3-D  configurations  for  later  predic¬ 
tions  of  aerodynamic  performance  and  to  the  subsequent  vis¬ 
cous  analysis  which  is  often  sensitive  to  minor  changes  in 
pressure  gradient.  The  result  is  that  the  proposed  method 
for  obtaining  uniformly  valid  solutions  about  arbitrary  con¬ 
figurations,  taking  into  account  the  effects  of  the  stagna- 
tion  regions,  may  be  a  useful  and  powerful  tool  in  the  aero- 
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dynamic  design  of  all  flight  vehicles. 


A  discussion  of  the  second  order  problem  is  presented 
in  Appendix  D.  A  necessary  and  sufficient  condition  is  pre¬ 
sented  for  the  existance  of  solutions  and  these  solutions 
are  shown  to  be  unique  to  within  an  unknown  constant.  The 
first  order  method  however  remains  preferable  to  a  second  or 
higher  order  effort  given  the  demonstrated  accuracy  of  the 
first  order  method  and  ease  of  application. 

Finally,  in  Appendix  E  interesting  comparisons  of  the 
term  as  obtained  for  the  2-D  ellipse  by  various 

approximations  and  the  exact  numerical  solution  are  shown. 
These  comparisons  further  demonstrate  the  validity  of  the 
first  order  method  in  obtaining  uniformly  valid  composite 
solutions  for  compressible  flow  about  bodies  with  round 
leading  or  trailing  edges. 
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II  The or 


The  purpose  of  this  chapter  is  to  provide  the  theoret¬ 
ical  background  necessary  before  general  application  of  the 
proposed  method.  The  exact  problem  will  be  stated  for  arbi¬ 
trary  configurations  in  potential  flow  with  a  discussion  of 
the  surface  boundary  c(  .xtion.  For  motivation  of  successive 
steps,  a  brief  order  of  magnitude  analysis  will  be  made  for 
the  exact  problem.  An  outer  expansion  will  then  be  assumed, 
identified  with  the  region  of  small  disturbances  from  the 
free  stream  conditions  for  the  flow  about  a  slender  body. 

The  resulting  first  and  second  order  approximations  to  the 
exact  potential  equation  will  then  be  stated.  Next,  an  in¬ 
ner  expansion  will  be  assumed  and  the  resulting  first  and 
second  order  approximations  for  a  stagnation  region  will  be 
identified.  Once  general  solutions  to  the  inner  and  outer 
problems  have  been  obtained,  one  would  classically  attempt 
to  match  these  solutions,  ultimately  forming  a  uniformly 
valid  composite  solution.  This  process  cannot  be  accomplish¬ 
ed  without  analytic  solutions.  Such  solutions  are  non-exis¬ 
tent  for  the  general  problem  (arbitrary  configuration)  leav¬ 
ing  as  one  alternative,  the  prospect  of  numerically  patching 
the  outer  and  inner  solutions  together.  To  avoid  this  pros¬ 
pect  the  theory  developed  at  this  point  is  reexamined  leading 
to  the  concept  of  an  intermediate  expansion.  This  expansion 
is  then  discussed  and  the  resulting  first  and  second  order 
approximations  are  identified.  With  the  intermediate  expan¬ 
sion  the  nature  of  a  uniformly  valid  composite  solution  is 
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hypothesized  and  then  formulated.  'i'he  result  is  a  method 
that,  can  be  applied  to  any  arbitrary  configure  tion  witn 
round  leading/trailing  edr.es  in  subsonic  potential  flow. 

The  first  order  method  as  developed  can  be  directly  applied 
using  any  numerical  technique  solving  the  Prandtl-Glauert 
equation.  This  versatility  leads  to  a  powerful  tool  for 
the  design  process,  especially  for  estimates  of  that  all 
important  performance  parameter,  drag. 

Statement  of  Exact  Problem 

The  steady  inviscid  irrotational  (potential)  flow  past 
an  arbitrary  configuration  (Fir  1 )  is  formulated  in  terms 


of  the  velocity  potential  $, 

where  the  velocity  V=y2> 

,  as 

follows  s 

\7M  )  --  -  [V 

•u(vdr] 

(11-1) 

pV-n  0 

on  surface  F(x,y,z)=0 

( I I -2 ) 

'  'J./C 

at  infinity 

( II  —  3 ) 

<-.*  «*  -  (■>  -  up 

(II-4) 

Figure  1. 


Axis  Syst.  t.  Notation 
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As  presented  here  the  surl'ace  boundary  condition  differs 

v 

form  l.ln'  u  l.niulii  rd  velocity  boundary  condition,  V-ri^O,  by  the 
inclusion  of  tho  local  donfsity.yo.  It.  is  reasonable  to 
assume  that  if  Eq  (ll-l)  could  be  solved  exactly  at  every 
point  in  the  flowfield,  then  since  the  density  is  not  zero 
anywhere,  it  would  suffice  to  enforce  the  standard  velocity 
boundary  condition  at  the  surface.  Equation  (11-1)  cannot 
be  solved  for  arbitrary  configurations  thus  one  must  resort 
to  an  approximate  method.  One  approach  is  to  attempt  using 
numerical  methods,  such  as  finite  difference  techniques,  to 
solve  Eq  (II-l)  with  the  appropriate  boundary  conditions  for 
a  given  configuration.  For  this  case  it  suffices  to  use  the 
standard  velocity  boundary  condition  at  the  surface  as  done 
by  Jameson  (Ref  16,17).  Another  approach  is  to  attempt  to 
solve  an  equation  which  approximates  Eq  (II-l),  The  pertur¬ 
bation  technique  employed  in  this  report  is  an  example  of 
this  latter  approach.  Equation  (II-l.)  is  a  mathematical 
statement  that  guarantees  conservation  of  mass  (continuity 
equation).  An  approximation  to  this  equation,  such  as  the 
brand t L-Glauert  equation,  does  not  assure  mass  conservation. 

It  has  been  claimed  by  Chin  (Ref  18)  that  to  obtain  analytic 
solutions  to  the  Prandtl-Olauert  equation  it  is  necessary  to 
reinforce  this  lost  continuity  condition  by  imposing  the  zero 
mass  flux  condition  at  the  surface.  In  fact,  it  appears  im¬ 
possible  to  transform  the  Prandtl-Glauert  equation  to  Laplace's 
equation  and  simultaneously  maintain  a  zero  net  flux  condition 
through  the  closed-body  streamline  (Ref  18)  when  trying  to 
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obtain  analytic  solutions,  unless  one  uses  a  consistent  ap¬ 
prox  Ilia  t  i  on  (  l.n  tho  m  pprox  i  m;i  t  i  .mi  implied  by  the  Fraud  LI. - 
Glauort  equation)  to  the  mass  flux  surface  boundary  condition. 
(This  is  demonstrated  in  Chapter  111  as  applied  to  the  elliptic 
cylinder. )  For  this  effort  then  the  mass  flux  surface  bound¬ 
ary  condition  and  approximations  thereto  will  be  used  through¬ 
out.  Further  demonstration  of  the  validity  of  using  the  mass 
flux  surface  boundary  condition  will  be  shown  in  Chapter  V. 
Finally,  Eq  (II-3),  is  the  standard  infinity  boundary  condi¬ 
tion  requiring  any  disturbances  due  to  the  presence  of  the 
body  to  disappear  at  fluid  infinity. 

The  exact  problem  can  be  restated  in  terms  of  a  pertur¬ 


bation  velocity  potential,  cp,  where 


$  -  Ur.  X  +-  (f  (II-5) 

The  gradient  of  the  perturbation  velocity  potential,  Vip,  thus 
represents  the  perturbation  velocity  components,  u,  v,  and  w, 
due  to  the  presence  of  the  body.  The  problem  can  now  be  re¬ 
stated  in  terms  of  non-dimensional  variables  as  follows* 

0  i  (.  x  <  p'  ^  ( 1 1 -6 ) 

<7>  -  Ul-7) 

X  =  y/i  ;  y  -  >'/(  ;  1  -  2/[i  ( II -S ) 

where  P  represents  some  characteristic 
body  length 

Replacing  the  local  speed  of  sound  in  Eq  (II-1)  with  Eq  (II-4) 
and  rewriting  Eqs  (II-l)  through  (1I-J)  in  terms  of  the  non- 
dimensional  variables  results  in  the  following  restatement 
of  the  exact  problem  (dropping  the  bar  notation  for  non-di¬ 


mensional  quantities) 
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Equation  (II-9)  dramatically  shows  the  highly  nonlinear  na¬ 
ture  of  this  problem.  Analytic  solutions  for  arbitraty  con¬ 
figurations  cannot  be  obtained  and  one  is  forced  to  seek 
numerical  or  approximate  solutions. 
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Order  of  Magnitude  Observations 

Before  attempting  to  develop  an  approximate  solution  to 
Eqs  (II-9)  through  (11-11),  some  interesting  observations  can 
be  made  by  a  rudimentary  order  of  magnitude  analysis  of  these 
exact  equations  (Ref  7).  The  classic  Prandtl-Glauert  equa¬ 
tion,  which  approximates  Ea  (II-9)»  results  to  first  order 
when  one  assumes  that  the  perturbations  due  to  the  body  are 
small : 


(■-A'  !W  A  ,/)  i- f/\  -  /s, 

\  '  1  a? '  f  a  x  v  *  *  A  ~  ( 


(11-12) 


Equation  (11-12)  forms  the  basis  for  thin  airfoil  and  slen¬ 
der  body  theories,  and  for  predictive  inviscid  aerodynamics 
today.  It  is  valid  for  slender  bodies  in  subsonic  or  super¬ 
sonic  potential  flows,  but  as  discussed  previously  yields 
incorrect  results  in  stagnation  regions.  If  one  examines 
Eq  (II-l)  or  (II-9)  for  the  stagnation  region,  where  V-0  and 
(or  -1),  to  first  order  Laplace's  equation  results: 

V  V  =  O  UI-13) 

Equation  (11-13)  describes  incompressible  potential  flow. 
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This  r«:;u  I  L  corves  ponds  wi  Lh  physical  Itiluilion  a:j  the  stag¬ 
nation  region,  however  small,  represents  a  region  whore  the 
total  velocity  and  Mach  number  are  approaching  zero,  thus 
approaching;  incompressible  (M=0)  flow  in  this  region  for  all 
free  stream  Mach  numbers.  Equation  (11-13)  while  yielding; 
the  correct  behavior  in  a  stagnation  region,  i3  clearly  in 
error  when  compressibility  effects  are  present  for  flow  over 
the  remainder  of  the  body.  Equations  (11-12)  and  (11-13) 
identify  the  distinctive  nature  of  the  flow  about  a  body  by 
establishing  two  different  regions  governed  by  two  different 
approximations  to  Eq  (11-9)  and  provide  the  motivation  for 
further  development. 

The  approach  taken  will  be  guided  by  singular  pertur¬ 
bation  theory  in  identifying  the  nature  of  these  regions  and 
establishing  the  appropriate  order  approximations  for  each. 

An  outer  expansion  will  be  assumed  associated  with  the  re¬ 
gion  of  small  disturbances  and  first  and  second  order  approx¬ 
imations  will  be  developed.  This  will  be  followed  by  the 
development  of  an  inner  expansion  associated  with  the  stag¬ 
nation  region.  The  analysis  will  be  restricted  to  round 
edges  in  subsonic  flow  as  it  is  felt  that  no  correction  is 
necessary  for  sharp  edges  (Ref  11,19).  The  restriction  to 
subsonic  flows  is  of  course  a  realistic  restriction  for  aero¬ 
dynamic  shapes  with  round  leading  edges.  The  method  as  de¬ 
veloped  could  also  be  applied  to  subsonic  round  leading  edges 
in  supersonic  flow. 

Outer  Expansion 


lo 


The  outer  ox  pans  ion  is  associated  with  tho  region  ol' 

Clow  about  n  body  where  the  aiuuimpl.i  on  of  :nmi  I  1  per  t.ui'bn  t,  i  onn 
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from  free  stream  conditions  is  valid.  This  is  of  course,  the 
basic  premise  of  slender  body  and  thin  airfoil  theory.  An 
outer  expansion  of  the  perturbation  velocity  potential  <p,  is 
assumed  in  terms  of  a  small  parameter,  £,  as  follows: 

< v*  -  £tftc  +  c-xq>2  +  -- "  (H-14) 

The  small  parameter,  G,  can  be  associated  with  either  the 
thickness,  camber,  or  angle  of  attack  of  the  body,  whichever 
is  appropriate  (Ref  20).  (For  later  applications  in  Chap¬ 
ters  III  through  V,  £  will  be  associated  with  thickness 
cases  only. )  Although  it  is  implied  by  the  form  of  Eq  (II- 
l4)  that  this  series  continues  indefinitely  in  powers  of  £, 
this  has  been  shown  to  be  not  necessarily  true  (Ref  19). 

For  a  round  nosed  profile  (2-D)  in  subsonic  flow,  logarithms 
of  £  beginning  with  6'V/:£  appear.  This  is  of  no  concern  for 
the  present  development  however,  as  terms  beyond  the  second 
order  will  not  be  used.  The  assumed  expansion  for  (fi°  (Eq  II- 
14)  incorporates  the  correct  behavior  to  be  expected  when 
there  is  no  disturbance,  or  as  C-K).  Substituting  Eq  (11-14) 
into  Eq  (II-9)  and  gathering  like  order  terms  with  respect 
to  the  parameter  C  results  in  the  following  first  and  second 


order  approximations  to  Eq  (11-9)  for  tho  outer  expansion: 

e  ■  ( '  -  M.i )  (,r  v y';  yy .« •  (^f/  =  o  ( 1 1  - 1 5 ) 

e}~  (i  -  -av*,  =  k!  U 

Equation  (11-15)  is  the  familar  Prandtl-Glauert  equation  as 


expected,  valid  for  subsonic  and  supersonic  flows. 
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A  surface  boundary  condi  Lion  is  certainly  appropriate 
for  this  expansion.  Assuming  the  isentropic  expression  for 
P/P-  and  small  perturbations  the  mass  flux  vector  can  be 

approximated  as  follows  (Ref  18,21): 

P  V  ~  /  i  ,  pi  R  v  w  \ 

f.  ...  '  ‘T/:..) 

The  surface  boundary  condition,  Eq  (II— 10),  then  becomes  (in 
terms  of  (^i0)  j 

(  I  *  .  *'  4'  *'■  »  “  O  on  P(x,y,z  )»0 

Substituting  Eq  (II-14)  into  the  preceding  equation  and  gafn 

ering  like  order  terms  of  the  parameter  6  results  in  the  fol 
lowing  approximations: 

I  :  '1  K  =  O 

(? '  /JVA'U  '/•;  -  >A"i  <r> 

6' :  //  A'  ^ .K  +  4Vy  A'  f  -  O 

The  zeroeth  order  approximation  (C°)  represents  a  flat  plate 
at  zero  angle  of  attack,  or  no  disturbances  (t=0),  as  ex¬ 
pected.  Now,  if  the  surface  is  described  by 

r(',y.A  -  cfiA  -  O 

where 

('  A'/  bd  . 

7  '  ,  fly  -  *'  /’• )  y/r- 1  >  0,  7  'vViVrl 

we  have  the  description  of  a  slender  body  at  zero  angle  of 
attack  with  e  representing  the  maximum  thickness.  The  sur¬ 
face  boundary  condition  now  becomes: 


"  i  <  dV,  '  ,'i  o 


'/  -y  O 


£  '  ■  aV;  ii  ■ 


4  'Ay  Sy  ‘  fj,  V4  - 


For  2-D  flows  the  first  order  approximation  yields 

■M  ,  . 

(p. ,  -  -  <>*/  •  "4  on  F(x,  z  )=0 
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For  configurations  with  round  leading  edges,  at  the 

n  lagan  t  i  on  point,  arid  t.lms  ^  t  would  bo  undofinod  in  this  re¬ 
gion  for  this  approximation.  It  is  proposed  then  that  the 
zeroeth  and  first  order  equations  of  the  original  approxima¬ 
tion  be  combined  as  follows: 

(  f  t  i-c-v/’/,  <•  on  F(x,  y ,  z  )=0  (11-17) 

This  then  forms  the  first  order  surface  boundary  condition 
and  is  similar  to  that  developed  by  Chin  (Ref  18).  The  sec¬ 
ond  order  surface  boundary  condition  approximation  is  then 
given  byt 

fi  VV.  c.  K  x  .'/V,  n,  -  o  on  F(x,y,z)-0  (II-18) 

Note  that  from  the  form  of  the  first  order  equation,  Eq  (II- 
17),  one  can  infer  immediately  that  the  stagnation  point 
velocity  will  not  be  zero  for  a  non-lifting  body  in  other 
than  incompressible  flow  (Mfr  =  0  or  /3  =1 ) .  This  again  confirms 
previous  statements  that  the  outer  solution  is  in  error  in 
this  region. 

The  infinity  boundary  condition  is  also  appropriate  for 
this  expansion  and  is  that  the  disturbances  due  to  the  pre¬ 
sence  of  the  body  disappear  to  first  and  second  order  at 
fluid  infinity: 

/yVy/_>r,  at  infinity 
/  ! -*  ^  at  infinity 

Inner  Expansion 

Before  assuming  an  inner  expansion  associated  with  the 
stagnation  region,  the  dependent  and  independent  variables 
previously  defined  are  strained  (Ref  10, iy)  to  magnify  this 


(11-19) 

(11-20) 
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region  where  the  outer  expansion  is  invalid.  Research  has 
shown  that  no  correction  is  noecssnry  ('or  sharp  edges,  cor¬ 
rections  will  bo  made  only  for  round  edges  where  the  associ¬ 
ated  radius  is  proportional  to  c*  (Ref  19).  Let  S-Cz  be  the 
parameter  used  to  strain  the  dependent  and  indepenoent  vari¬ 
ables  as  follows: 

(]>  --  */<,  •  x  -  */s  ;  y  --  tk  ■  Z/Z  (11-21 ) 

Equation  (II-9)  remains  unchanged  in  terms  of  the  new  vari¬ 
ables.  An  inner  expansion,  (p'd.y,?.)*  is  then  assumed  in  powers 
of  where, 

(() 1  -■  -  x  +  f  (  ••  *  -  (11-22) 

The  reason  for  this  choice  becomes  clear  as  one  sees  that  as 
the  disturbance  parameter,  S,  approaches  zero  (flat  plate  at 
zero  angle  of  attack),  the  total  velocity  approaches  zero  or 
stagnation  point  conditions  as  desired.  Substituting  Eq  (II- 
22)  into  Eq  (II-9)  in  terms  of  the  new  variables  results  in 
the  following  first  and  second  order  approximations: 

(11-23) 

S7'  \) •vy. ;  :0  (11-24) 

(The  Laplacian  operator  denoting  operations  with  respect  to 
x,  y,  and  z. )  Rewriting  Eqs  (11-23)  and  (11-24)  in  terms  of 
the  original  variables  results  in  the  following: 

t  :  \J  V*;  -  (")  (11-29) 

<\  \7 -  O  (II-2o) 

(The  Laplacian  operator  now  denoting  operations  with  respect 
to  x,  y,  and  z. )  Equations  (11-25)  and  (II-2b)  again  con¬ 
firm  physical  intuition  and  the  order  of  magnitude  analysis 
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presented  earlier  that  the  1‘low  in  the  stagnation  region  can 
bo  approximated  as  being  incompressible,  even  to  second  order 
(this  is  not  true  for  higher  order  approximations).  This 
implies  that  a  solution  can  be  obtained  independent  of  free 
stream  conditions  once  the  shape  is  specified,  for  subsonic 
through  supersonic  flows  (assuming  the  appropriate  boundary 
conditions  have  been  specified). 

A  surface  boundary  condition  is  also  appropriate  for 
this  region.  To  second  order,  an  expansion  of  p/p ^  from 
stagnation  conditions  is  a  constant  as  expected.  Since  the 
flow  in  this  region  is  being  approximated  by  incompressible 
flow  (M=0),  the  surface  boundary  condition  is  identical  to 
the  velocity  boundary  condition  at  the  surface  (written  in 
outer  variables )  : 


( 1 1  ;r 


..l  \ 


,'y.  +  --n 


on  F(x,y,z)=0 


Substituting  Eq  (11-22)  into  the  preceding  equation  in  terms 
of  inner  variables  and  gathering  like  order  terms  of  the 
parameter  &  results  in  the  following  first  and  second  order 
surface  boundary  condition  approximations  (written  in  outer 
variables ) : 

v/V\  -  *1  -  C>  on  F(x,y,z)=0  (11-27) 

\]p\  -■  /  -  O  on  F(x,y,z)  =  0  (11-28) 


Summary 

As  noted  earlier  the  outer  solution,  once  obtained,  is 
valid  everywhere  about  a  slender  body  or  thin  airfoil  except 
in  the  stagnation  region.  An  inner  solution  would  of  course 
be  valid  only  in  the  stagnation  region.  It  appears  attractive 
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and  desirable  to  combine  these  two  solutions,  retaining  the 
OGUonllal  features  of  each,  forming  a  uniformly  valid  com¬ 
posite  solution.  Classically  this  would  be  done  (to  first 
order  for  example)  by  first  solving  the  outer  problem  given 
by  Eq  (11-15)  with  boundary  conditions  given  by  Eqs  (11-17) 
and  (11-19)  s  and,  the  inner  problem  given  by  Eq  (11-25)  with 
the  boundary  condition  given  by  Eq  (II-2?)  to  within  an  un¬ 
known  function.  Next,  the  solutions  would  be  matched  (Ref 
19)  to  determine  the  remaining  boundary  condition  and  form 
a  uniformly  valid  composite  solution.  This  composite  solu¬ 
tion  is  obtained  by  summing  the  outer  and  inner  problems  or 
solutions  to  appropriate  order  and  subtracting  their  common 
limit.  This  process  requires  analytic  solutions  and  as  men¬ 
tioned  previously,  such  solutions  are  nonexistent  for  appli¬ 
cations  to  arbitrary  configurations.  The  remaing  alternative 
is  to  numerically  patch  the  two  solutions  together.  As  the 
first  approach  is  not  plausible,  this  approach  is  not  desir¬ 
able.  The  problem  remains:  how  to  form  a  uniformly  valid 
solution  in  a  reasonable  manner.  Before  accomplishing  this 
however,  the  ability  to  match  outer  and  inner  solutions  as¬ 
sumes  that  their  respective  domains  of  validity  overlap. 

There  is  no  reason  a  priori  to  assume  that  this  is  true  and 
in  fact  for  some  applications  they  may  not  overlap  (Ref  19). 

Intermediate  Expansion 

To  bridge  the  gap  that  may  exist,  between  outer  and  inner 
domain:;  of  validity  it  is  necessary  to  consider  an  interme¬ 
diate  expansion,  a  concept  developed  mainly  by  Kaplun  (Ref 
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19.22).  It  was  introduced  in  part  to  try  and  formally  ex¬ 
plain  why  ihe  classical  method  of  matching  solutions  worked 
in  cases  where  it  should  not,  and  did  not  work  in  cases 
where  it  should.  Kaplun's  Extension  Theorem  (Ref  19,22) 
attempts  to  formally  define  the  matching  process  through  de¬ 
fining  and  extending,  if  necessary,  regions  or  domains  of 
validity  by  using  the  intermediate  expansion.  This  attempt 
is  viewed  as  heuristic  by  some  (Ref  19).  To  begin,  the  de¬ 
pendent  and  independent  variables  are  again  strained  as  fol¬ 
lows  t 


A/  ( 0/ 

V  -  'Y-:a  ;  x  ■-  /<?'•;  y 


%*  :y  -  V(-A  ;  1  --  */e 


(11-29) 


Equation  (II-9)  remains  unchanged  in  terms  of  the  new  vari¬ 
ables.  Letting  X  -€C%  an  intermediate  expansion,  Cpl(pi  y  i\  is 
assumed  as  follows* 


cf1  --  -  ?;  X  1  r,j V  «  t'upf 


A 


(11-30) 


The  reason  for  this  choice  can  be  clearly  seen  since  for 
A=0  the  outer  variables  result,  while  for  A=2  the  inner  vari¬ 
ables  and  expansion  result.  Thus  A  is  assumed  to  be  some  un¬ 
determined  constant  between  zero  and  two,  i.e., 


0  <  A  <  2  (II-31) 

This  expansion  can  be  physically  interpretted  as  an  expansion 
from  while  the  outer  expansion  is  from  U,u ,  and 

the  inner  expansion  is  from  a  zero  stagnation  point  velocity. 
Substituting  Eq  (II-30)  into  Eq  (II-9)  in  terms  of  the  new 
variables  results  in  the  following  first  and  second  order 
approximations  (rewritten  in  the  original  variables): 

*  -.A.*}"  fp.l  .1  r  O  (11-32) 


. .a* • .  afc Aa 
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Again  note  that  for  A~0  t.ho  outer  equation:*,  result  wnilo  for 

A=2  the  inner  equation'*,  result. 

A  surface  boundary  condition  also  seems  appropriate  for 
this  region.  Physically,  the  first  order  intermediate  approx¬ 
imation,  Eq  (11-32),  appears  to  be  flow  at  a  Mach  number  de¬ 
fined  by  Mt -H-A,  since  A  appears  in  the  same  position  as  /3J 
in  Eq  ( 1 1  —  1  *3 )  for  the  outer  approximation  (A  yV’  for  A  =  0). 

Based  on  this  analogy  the  surface  boundary  condition  for  the 
intermediate  approximation  is  assumed  to  be: 

(  I  i  AtpA/l.  +  <l\LA,  4-  --  o  on  F(x,y,  z  )  =  0 

In  terms  of  the  outer  variables  the  first  and  second  order 

boundary  condition  approximations  then  are: 

[  A  -  §  A')  4-  tht>\  j -1, 4  -1,  e  t->)  %-0  on  F  ( x ,  y ,  z  )  =  0  (11-34) 

■Af/VyA  -  O  on  F(x,y,z)  =  0  (11-35) 

Note,  as  before,  that  for  A=0  the  outer  equations  result 
while  for  A=2  the  inner  equations  result.  Now  by  using  the 
intermediate  expansion  and  applying  Kaplun's  Extension  Theo¬ 
rem  (Ref  22),  the  inner  and  outer  solutions  may  be  formally 
matched.  This  is  accomplished  essentially  by  extending  the 
domains  of  validity  of  each  (if  necessary)  and  choosing  A 
such  that  these  domains  overlap.  As  before,  however,  this 
is  a  process  which  appears  to  require  analytic  solutions. 


Hypotheses 

The  problem  remains  as  to  how  to  form  a  uniformly  valid 
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composite  solution.  To  avoid  numerical  patching  it  is  neces- 

Gary  to  somehow  detorinme  the  common  limit  required  to  form 

this  composite  solution.  Consider  the  following  quotes  with 

respect  to  the  composite  solution  and  intermediate  expansion* 

"Composite  solution. .. the  sum  of  the  inner  and  outer 
expansions  is  corrected  by  the  part  they  have  in  common 
(Ref  19*94). " 

"Intermediate  expansion. .. in  some  overlap  domain  the 
intermediate  expansion  of  the  difference  between  the 
outer  (or  inner)  expansion  must  vanish  to  appropriate 
order  (Ref  19*92)." 

In  light  of  the  essence  of  these  statements  and  the  inherent 
nature  of  the  intermediate  expansion  with  respect  to  the 
outer  or  inner  expansions,  the  following  hypothesis  is  pro¬ 
posed  * 

i.  The  solution  to  the  intermediate  expansion  is,  to 
sufficient  order,  the  common  limit  necessary  for 
^  matching. 

With  this  hypothesis,  a  composite  solution  can  now  be  formed. 
However,  the  inner  and  intermediate  problems  cannot  as  yet 
be  fully  solved  as  they  lack  a  necessary  boundary  condition. 
Furthermore,  for  application  to  current  numerical  methods, 
it  would  be  desirable  to  have  each  problem  (outer,  inner,  and 
intermediate)  be  independent,  well  posed  and  unique  to  a  giv¬ 
en  configuration.  Each  solution  contributing  to  the  forma¬ 
tion  of  a  uniformly  valid  composite  solution  as  to  be  demon¬ 
strated.  Since  the  nature  of  thiG  problem  is  a  uniform  stream 
perturbed  by  a  body  whose  disturbances  must  disappear  at  in¬ 
finity,  a  final  hypothesis  is  proposed* 

ii.  The  remaining  boundary  condition  for  the  inner  and 
— ^  intermediate  problems  is  that  disturbances  due  to 

'  the  body  disappear  at  fluid  infinity. 
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■v  J  t  h  hy  do  these:;  (i)  and  (ii>  t.hree  well  posed  problems  are 
defined,  each  ol"  whicn  can  bo  solved  irule penden lly  for  the 
same  configuration,  such  that  each  solution  contributes  to 
a  uniformly  valid  composite  solution. 


r> 


Composite  Solution 

As  previously  stated,  the  composite  solution  is  the  sum 
of  the  outer  and  inner  solutions  minus  the  part  they  have  in 
common.  As  hypothesized,  this  common  limit  is  the  interme¬ 
diate  solution  for  a  given  value  of  A.  The  resultant  veloc¬ 
ities  for  each  solution,  being  directly  a  function  of  the 
individual  perturbation  potential  functions,  will  then  be 
calculated  and  used  to  form  the  composite  solution  (Ref  19) • 
The  use  of  the  respective  velocities  to  form  the  composite 
solution  is  a  matter  of  choice  as  a  composite  could  just  as 
easily  be  formed  in  terms  of  the  perturbation  potential  <p. 
The  individual  velocities  however,  are  more  definitive  in 
terms  of  physical  significance  as  well  as  being  direct  in¬ 
puts  to  any  aerodynamic  calculations.  The  velocities  for 
each  region  in  terms  of  the  perturbation  potentials  for  the 
respective  expansions  are  given  by  the  following: 

^  '  (  1  .  Vv  -  <i'‘t  ) 

-0  *  *  fvs ,  e.4'',  * 


V' 


'/■;  -  /; ) 


u«.  ~  i  "4' 

-  ( '  'i  i*  <  V.V  *  -s'vVr ,  -V/.j  r  yrj. ) 

i  -■ 

Z(  I >  trf >  t ;  ,  ti,  :‘,i  V  v ) 


(11-30 ) 


(11-37) 


(11-33) 


The  composite  solution  may  be  formed  by  either  of  two  combi- 
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nations  (Kef  1  V )  oi'  the  resultant  velocities  from  each  re- 
... 

gion  as  follows  (nole  that  lor  incompressible  How,  ^=0, 

both  rules  yield  the  same  result  as  desired  )i 

VI  Vi  vi  vl 

(additive)  ],\f  ‘  :j.Jc  '  u, !,  "  i/..  I  i  Ul-39) 

(multiplicative)  '^|f  j  (  \})[  )}/'  V  j  ^  (lI-40) 

In  terms  of  the  surface  velocity,  the  multiplicative  compos¬ 
ite  solution  is  preferred,  as  it  has  been  shown  (Ref  19*23) 
to  give  better  results  in  the  stagnation  region  (this  is 
easily  seen  by  the  fact  that  the  multiplicative  composition 
results  in  a  zero  velocity  at  the  stagnation  point  for  all 
free  stream  Niach  numbers).  Note  however,  that  in  using  the 
multiplicative  composition,  the  possibility  exists  for  the 
intermediate  solution  to  equal  zero,  thus  resulting  in  a 
singular  point.  This  has  been  shown  (ref  23)  to  occur  in 
the  flowfield  and  not  on  the  surface,  and  will  be  verified 
later  analytically. 

Summary 

The  following  is  a  summary  of  the  first  and  second  or¬ 
der  approximations  for  each  expansion  as  developed  in  this 
chapter. 

i.  Outer  Expansion: 

<t':„  -o  (ii-4i) 

(i 1  f  df/yl,,  Mu/pd,  0  on  F(x,y,z)=0  (11-42) 

W,/~>  O  at  infinity  (11-43) 

^V/Vk.tc/i;v bail,  l 

+  2  til  jf/Vr  L<]ly  1  </?>  <i:u  J  ( 1 1  -44  ) 

^  'I  *- !/[-.  Vy  t  --  O  on  F(x,y,z)=0  (11-45) 

I  (/ tfl  I  -  ■>  ("■  at  infinity  (II  -4b  ) 
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ii.  inner  Expansion: 

(11-4?) 
(11-48) 

(!'[. -b h.  j-'O  n  1  J  ni  lnny  (11-49) 

(11-50) 
(11-51) 
(11-52) 

Intermediate  Expansion: 

ML  --O  (11-53) 

[(/-  MA-0  on  F(x , y , z )=0  (11-54) 

(f.[  ~'>  >  A'!  Aw! )  ->  C  at  infinity  (11-55) 

.AcA;,  *  c/ly,  4  '/Av  -in  !  t 

t  -  A,  (11-56) 
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A(A:  /?« 95:  A  *  vlv  A  =  O  on  F(x,y, z )=0  (II -5?  ) 

/V'Av/-^0  at  infinity  (II-58) 

The  composite  solution,  in  terms  of  resultant  velocities,  is 
formed  according  to  Eq  (11-39)  or  (11-40).  It  is  interest¬ 
ing  to  note  that  all  the  first  order  problems  as  stated  can 
be  solved  using  current  numerical  techniques,  such  as  panel 
methods  (Ref  5)»  which  solve  the  Prandtl-G lauert  equation 
subject  to  mass  flux  surface  boundary  conditions.  The  outer 
problem  would  be  solved  for  the  free  stream  Mach  number,  the 
inner  problem  for  a  zero  Mach  number  (incompressible),  and 
i;he  in  termed  into  problem  for  some  intermediate  Mach  number, 

Mj,  as  determined  by  the  value  of  A.  The  more  exact  and  uni¬ 
formly  valid  composite  solution  would  then  be  formed  as  speci¬ 
fied  making  the  first  order  method  a  very  powerful  yet  simple 
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technique  1'or  application  to  arbitrary  configurations  in  sub¬ 
sonic  potential  flow. 
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The  purpose  of  this  chapter  is  to  describe  the  detailed 
results  obtained  by  applying  the  first  order  approximation, 
as  developed  in  Chapter  II,  to  a  nonlifting  elliptic  cylin¬ 
der  (2-b  ellipse)  in  subsonic  compressible  flow.  The  ellip¬ 
tic  cylinder  was  chosen  for  application  so  that  analytic 
rather  than  numerical  results  coulc  be  used  to  test  the  meth¬ 
od  for  a  variety  of  thicknesses  and  Mach  numbers.  The  pro¬ 
blem  will  be  stated  with  the  surface  mass  flux  boundary  con¬ 
dition  and  an  analytic  solution  will  be  presented.  Details 
of  the  solution  are  in  Appendix  A.  This  is  the  first  such 
solution,  to  this  author's  knowledge,  with  the  prescribed 
boundary  conditions.  The  necessity  for  the  mass  flux  bound¬ 
ary  condition  will  also  be  demonstrated  as  promised  earlier. 
Surface  velocities  will  be  evaluated  for  an  ellipse  with  a 
ten  percent  maximum  thickness  to  chord  ratio  at  a  free  stream 
Mach  number  of  0.5.  This  case  will  be  used  to  depict  the 
various  aspects  of  the  proposed  method  and  solutions,  and 
show  that  the  multiplicative  composition  is  preferred  to  the 
additive  composition  for  representation  of  surface  velocities. 
No  attempt  will  be  made  to  optimize  the  value  of  A ,  Eq  (II— 
31),  and  unless  otherwise  stated,  the  value  of  A  will  be  one 
for  all  applications.  Finally,  the  first  order  multiplica¬ 
tive  composite  results  will  be  compared  to  results  from  thin 
airfoil  theory  (Appendix  B)  and  numerical  results  obtained 
from  a  Jameson  code  (Ref  16,17)  which  solves  the  exact  poten- 


Li;il  equation  i  L'.'q  (.11  — v  J  with  the  veioei  Ly  suri'aco  boundury 
condition  and  l, ho  standard  infinity  boundary  condition). 
Comparisons  will  bo  made  of  total  surface  velocity  and  sur¬ 
face  velocity  perturbation  components  (u  and  w)  for  a  vari¬ 
ety  of  thicknesses  and  Mach  numbers.  'lhe  comparisons  will 
show  that  the  uniformly  valid  first  order  solutions  compare 
very  favorably  with  the  exact  numerical  solutions  within  the 
limits  of  validity  of  the  theory,  that  is,  for  thin  bodies 
at  subsonic  Mach  numbers.  The  composite  results  do  not  com¬ 
pare  favorably  at  transonic  Mach  numbers  or  for  thick  bodies 
in  subsonic  or  transonic  flow. 

Statement  of  Problem 

For  application  of  the  first  order  method,  consider  a 
nonlifting  elliptic  cylinder  (Fig  2)  in  the  subsonic  com¬ 
pressible  isentropic  flow  (potential  flow)  of  a  perfect  gas 
(Y=1.4),  whose  surface  is  given  by 

P (*./.)  =  xz  i-  -J  =  0  (III-l) 

Here  the  parameter  £  is  the  maximum  thickness  to  chord  ratio. 


Figure  2.  Flliptic  Cylinder 


The  outer,  inner,  ami  intermediate  first  order  problems  are 
all  of  the  form 


i Yq\^  icA„  -O  (III-2J 

f  on  ?(x,z)=0  ( 1 1 1 -3  ) 

(f:  K  >  ^  r  ;  </'■,  j  — »  O  at  infinity  (III-4) 

where 

fo2,  outer  M  »  2"  t<t 

tj“- |A  ,  urtterejgfJidte  ( <3< 4  <-2  ),  c  .•  i- f  ^  (III-5) 

(J  ,  mi  her  ((\:  2) ,  Z  -A1 

Equation  (I1I-2)  with  boundary  conditions  given  by  Eqs  (III- 
3)  and  (III-4)  states  the  first  order  problem  for  an  elliptic 
cylinder.  The  specific  approximation  is  determined  by  the 
value  of  A,  Eq  (II1-5). 

Method  of  Solution 

Z"''  The  solution  method  (details  in  Appendix  A)  for  obtain¬ 

ing  analytic  results  for  the  first  order  problem  requires 
transformation  of  the  problem  to  a  known  solution,  incompress¬ 
ible  flow  about  a  circular  cylinder.  This  involves  first 
transforming  the  problem  to  incompressible  flow  about  an 
equivalent  elliptic  cylinder  with  the  transformations 
and  =  ('-  irY)\  +  qi  resulting  in: 

uu-6) 

(TxVf  3  -C.  on  K(x,z)=0  (III-?) 

•  )  ‘Lit-’  at  infinity  (111-6) 

where 

hd.i)  -  L'BxY't  (y'Az  -  /  =c>  (m-9) 

Once  this  is  accomplished,  the  problem  is  then  mapped  using 

/T' 

elliptic  coordinates  to  the  known  solution.  The  mapping 
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function  used  is  t.he  well  known  joukowski  Iraiui l'orm, 


The 


/T 


resulting  first  order  analytic  solution  is  then, 

vrtu . *\  •  !f$  ( (? * k'- 1’ ] [ i  *  w'} -  (i-  l]f  Vo'  '  °  n'J- hU 


where 


k  =  ~  fT^TKF 

Q  --  VV4  k 

p  -  *1  i  ( Yp  *  kY 
<i  ■-  fik  *• ;  vo  kp 


and 


(111-10 j 

(III-U) 
(III  — 12 ) 
(III -13) 
(III-14) 
( III -1 5 ) 
(III-16) 


u/tL.  -cA,  ,  =  (III-17 ) 

r  (^y'  (III-18) 

For  a  given  free  stream  Mach  number,  MCi  ,  and  maximum  thick¬ 
ness  to  chord  ratio,  €,  the  velocities  for  each  solution  can 
be  evaluated  with  A=0  for  the  outer  solution,  A=2  for  the 
inner  solution,  and  an  assumed  value  of  A  (A=l  )  for  the  inter¬ 
mediate  solution.  The  velocities  for  the  composite  solution 
are  then  determined  using  either  Eq  (11-39)  or  (11-40). 

Let  us  now  determine  what  would  have  happened  if  the 
standard  velocity  surface  boundary  condition,  and  not  Eq 
(III -3 ) .  had  been  applied.  Consider  the  following  restate¬ 
ment  of  the  problem  with  the  velocity  boundary  condition, 

V  n-Q,  being  applied  at  the  surface.  The  infinity  boundary 
condition  is  not  affected. 

0  V...  '  <t\,,  '  V. 

O-J  i  i-  (m-19) 
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(Note  that  a  xero  s  lugna  Lion  point  velocity  is  enforced  as 
a  consequence  of  using  Kq  (111-1V)  for  all  free  stream  Mach 
numbers  for  the  outer  approximation,  A  =  0.  We  must  ask  if 
this  is  a  reasonable  result  of  a  method  based  on  small  per¬ 
turbations,  outer  approximation,  for  compressible  flow.) 

The  standard  Prandtl-Glauert  transf ormation,  x=Bx,  is  applied 
resulting  in, 

<P. «  *  '/'}«  -C-> 

(I  ~  '}  4  t-  '-(/.  Ac1)  r  O  (III-20) 

Now  let  (p and  (p~- *■  'J!  resulting  in 

I'n;  '  <?*,*,  -O 

GW  cl.AA)  --O  (HI-21) 

As  before,  the  problem  has  been  successfully  transformed  to 
the  incompressible  flow  case,  but  now  the  surface  boundary 
condition  is  being  satisfied  on 

t\  --  xz  AA-  Y  -  l  =  C-  (Ill -22 ) 

or 

/-  (x. it)  -  (  ApA  1-  ( -  !  =  O  ( III-23 ) 

and  not  on  the  specified  surface,  F(x,z)=0,  given  by  Sq 
(III-l).  Thus  while  one  may  obtain  a  solution  that  appears 
correct,  since  both  surfaces  have  the  same  maximum  thickness 
(only  the  chord  lengths  are  different)  and  yield  approxi¬ 
mately  tho  same  maximum  velocity,  it  is  clearly  for  the  wrong 
surface.  In  fact,  this  problem  cannot  be  successfully  trans¬ 
formed  to  an  equivalent  incompressible  flow  problem  for  the 
same  or  given  surface  using  the  velocity  boundary  condition 
(Ref  18).  Recall  that  thin  airfoil  theory  neglects  the  x- 
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component,  of  perturbation  velocity  and  thuo  this  problem  does 
not  arise.  This  results  for  any  surface  and  is  not  peculiar 
to  the  ellipse. 

Analytic  Results 

Analytic  results  of  the  first,  order  method  will  now  be 
presented  and  discussed  for  the  2-D  ellipse.  The  results  are 
for  an  ellipse  with  a  ten  percent  maximum  thickness  to  chord 
ratio  at  a  free  stream  Mach  number  of  0.5.  All  data  is  for 
the  nonlifting  problem  only. 

Figure  3  shows  the  surface  velocity  comparisons  for  the 
outer,  inner,  and  intermediate  (A=l )  solutions.  These  three 
solutions  of  course  contribute  to  the  composite  solution. 

In  Fig  3  it  is  clearly  seen  that  the  outer  solution  is  in 
error  at  the  stagnation  point.  The  inner  solution,  although 
correct  at  the  stagnation  point,  does  not  include  the  effects 
of  compressibility  as  represented  by  the  outer  solution  for 
the  remainder  of  the  body.  The  intermediate  solution  is  seen 
to  be  in  error  for  both  of  these  regions,  N'one  of  these  in¬ 
dividual  solutions  represents  a  solution  that  can  be  consid¬ 
ered  uniformly  valid  for  the  entire  body. 

Figure  4  presents  the  surface  velocity  comparisons  for 
the  outer,  inner,  and  multiplicative  composite  (A=l)  solutions. 
Here  it  is  seen  that  the  composite  solution  combines  the 
compressibility  of  the  outer  solution  with  a  zero  stagnation 
point  velocity  representative  of  the  inner  solution.  The 
composite  solution  transitions  smoothly  from  the  inner  to 
the  outer  solution,  combining  the  effects  of  both.  The  tran- 


sillon  occurs  very  ciojo  to  the  stagnation  point,  confirming: 
the  fact  that  the  outer  solution,  or  Prand tl-Glauert  solu¬ 
tion,  is  the  dominant  solution  (and  as  previously  discussed, 
a  very  good  approximation,  borne  out  by  its  common  use,  ex¬ 
cept  in  the  stagnation  region).  The  pertinent  result  in 
this  figure,  however,  is  that  the  composite  solution  repre¬ 
sents  a  uniformly  valid  solution  about  a  thin  body  with  blunt 
edges. 

Figure  5  shows  the  effect  when  using  the  additive  com¬ 
posite  solution  on  the  surface  velocity  due  to  varying  the 
value  of  A  (A=0. 5, 1. 0, l. 5  ) .  As  A  increases  more  of  the  ef¬ 
fect  of  the  inner  solution  is  being  subtracted  from  the  com¬ 
posite  solution,  and  the  composite  solution  moves  toward  the 
outer  solution  as  shown.  The  difference  between  the  solu¬ 
tions  for  the  different  values  of  A  is  negligible  except  at 
and  near  the  stagnation  point.  The  additive  composition  was 
used  here  to  illustrate  these  effects  and  also  to  show  why 
the  multiplicative  composition  is  preferred  since  a  zero 
stagnation  point  velocity  results  from  the  multiplicative 
composition  for  all  values  of  A.  This  comparison  between 
the  additive  and  multiplicative  compositions  is  shown  in  Fig 
6  for  a  value  of  A=l.  As  shown  in  this  figure  the  solutions 
are  identical  except  at  and  near  the  stagnation  point.  These 
figures  suggest  that  there  may  be  some  optimum  or  desired 
value  of  A  for  forming  the  composite  solution,  however,  this 
requires  some  other  criteria  than  has  been  presented  (such 
as  experimental  data).  Without  such  guidance,  it  appears 
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that  a  value  of  A;  l  ilioultl  uui'fLoo  for  mot;  t.  application:;. 

A::  mon  l  i  on«M  n:irl  iim  ,  Un>  puraii  P  i  l  i  l.y  exist:;,  whoti  using 
t.hn  iriu  1  t,  i  p  l  i  ca  (.  i  vo  composition,  that,  (ho  i  n  to  rrnnd  i  n  to  aolu- 
tion  may  equal  zero  at  some  point  yielding  a  singular  solu¬ 
tion  at  this  point.  This  does  indeed  occur,  as  shown  in  Fig 
7,  but  always  occurs  in  the  flowfield  and  not  on  the  surface 
(Ref  23).  This  demonstrates  that  although  the  multiplicative 
composition  is  preferred  for  evaluating  surface  velocities, 
care  must  be  taken  when  using  this  composition  in  the  flow- 
field.  The  additive  composition  does  not  have  any  singular 
points . 

The  purpose  of  tne  figures  presented  to  this  point  has 
been  to  show  that  the  three  first  order  analytic  solutions 
(outer,  inner,  intermediate)  can  be  successfully  combined 
into  a  uniformly  valid  solution.  This  composite  solution 
correctly  combines  compressibility  effects  (represented  by 
the  outer  solution)  with  the  correct  flow  behavior  in  the 
stagnation  region  (represented  by  the  inner  solution)  for  a 
thin  body  with  blunt  edges  in  subsonic  flow.  The  composite 
solution  is  formed  by  combining  these  two  solutions  being 
sure  to  subtract  their  common  limit  (represented  by  the  in¬ 
termediate  solution  for  a  fixed  value  of  A).  The  three  solu¬ 
tions  are  essentially  solutions  to  the  Prandtl-Glauert  equa¬ 
tion  for  (outer  solution),  and  r."=Mx  (intermediate  solu¬ 

tion  where  M(<<Mt<0),  and  M=G  (inner  solution),  for  the  given 
surface  and  infinity  boundary  conditions.  In  terms  of  sur¬ 
face  velocities  the  multiplicative  composition  has  been  shown 
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SURFACE  VELOCITY  (v/U 


Figure  3.  First  Order  Surface  Velocity  Comparisons, 
Outer,  Inner,  and  Intermediate  (A=l  )  Analytic  Solutions, 
2-D  Ellipse,  €=0.1,  Mco=0.'5,  0-50;i  Chord 


Figure  4.  First  Order  Surface  Velocity  Comparisons, 
Outer,  Inner,  and  Multiplicative  Composite  (A  =  l )  Analytic 
Solutions,  2-D  Ellipse,  €=0.1,  M^O.5,  0-50$  Chord 
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Figure  5.  First  Order  Stagnation  Point  and  Surface 
Velocity  Comparisons,  Outer,  Inner,  and  Additive  Compo: 
ite  ( A=0.  5, 1 . 0, 1. .  << )  Analytic  Solutions,  2-D  Ellipse, 
£=0.1,  1^=0.  5,  0-50#  Chord 
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Figure  6.  First  Order  Surface  Velocity  Comparisons, 
Additive  and  Multiplicative  Composite  Analytic  Solutions 
^  (A  =  l ) ,  2-D  Ellipse,  C=0.1,  M^O.  S,  0-50#  Chord 
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VELOCITY  (V/U, 


Figure  7.  Singular  Point  Location,  First  Order  Analytic 
Multiplicative  Composite  Solution  (A=l),  2-D  Ellipse,  £=0.1, 
^=0,  5 
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to  be  preferred  to  the  additive  composition.  Next,  the  first 
order  nm  I  I.  i  p  I  i  ca  I  i  v<»  composite  i’i'mhI  I.::  will  bo  compare*.!  to 
the  results  obtained  from  other  methods  for  a  variety  of  free 
stream  Mach  numbers  and  thicknesses  for  the  ellipse. 

Comparison  with  Other  Methods 

The  best  comparison  that  could  be  made  would  be  with  an 
exact  analytic  solution  to  the  full  potential  equation,  Eq 
(II-9).  This  is  not  possible  and  so  the  results  are  compared 
to  numerical  solutions  obtained  from  a  Jameson  code  (Ref  16, 

17)  modeling  the  full  potential  equation.  The  first  order 
multiplicative  composite  results  are  also  compared  with  the 
classic  analytic  results  from  thin  airfoil  theory.  The  com¬ 
parisons  will  be  made  in  terms  of  the  total  surface  velocity, 
V/U m  (first  order  composite,  Jameson,  and  thin  airfoil  theory), 
and  surface  perturbation  velocity  components,  u/u^  and  w/u^ 
(first  order  composite,  Jameson).  The  comparison  will  be 
for  ellipses  of  £=0.8,  at  Ma.=0,l,  £=0.5,  at  M^O.  5,  and  £=0.1, 
at  M^O. 1,0. 3, 0.5,0. b,0.7,  and  0.8,  It  should  be  noted  that 
the  results  from  thin  airfoil  theory  represent  an  asymptotic 
value  in  terms  of  the  maximum  surface  velocity  for  the  first 
order  method  (this  is  consistent  with  the  derivation  and  the 
theory ) . 

Figures  8a  through  8f  (legend  for  Figs  8a  through  8f  in 
Table  I)  present  the  surface  velocity  and  perturbation  velo¬ 
city  components  for  the  first  order  multiplicative  composi¬ 
tion  (A=l),  exact  numerical,  and  thin  airfoil  theory  solu¬ 
tions.  This  data  is  for  an  elliptic  cylinder,  £=0.l,  at  free 
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Mach  number:!  of  0.1,0.  ) ,  O.  u.  o,  0. /,  and  0.13,  The 
(lain  at,  M<0  o.  1  are  oar.  on  t.  La  !  I  y  i  neompress  i  b  Lo  while  the  data 
at  M^O.8  are  near  the  transonic  regime.  (For  most  applica¬ 
tions  M^O.3  is  the  limit  for  approximating  the  flow  as  in¬ 
compressible  while  the  transonic  flow  regime  is  approximately 
0. 8  <  Mct<  1 . 2.  'these  of  course  are  only  guidelines  and  depend 
on  the  particular  configuration. )  The  first  order  multipli¬ 
cative  composite  solutions  are  seen  to  compare  very  favorably 
with  exact  numerical  solutions  at  the  lower  Mach  numbers  ( Mw= 
0.1, 0.3,  and  0.5).  As  Mach  number  increases,  the  agreement, 
although  still  very  good,  begins  to  decrease  until  at  M^O.8 
the  composite  results  are  seen  to  be  clearly  in  error.  The 
discrepancies  are  due  to  the  dominance  of  the  outer  solution 
in  the  composite  solution.  For  the  maximum  velocity,  this 
is  primarily  due  to  the  fact  that  as  Mach  number  increases 
from  subsonic  to  transonic  Mach  numbers,  the  assumption  that 
the  perturbations  remain  small  becomes  less  valid.  This  is 
clearly  shown  by  the  thin  airfoil  theory  solution,  being  the 
asymptotic  value  of  the  first  order  composite  result.  This 
disagreement  between  the  composite  and  exact  solutions,  as 
Mach  number  increases,  seems  to  arise  primarily  in  the  x- 
component  of  the  perturbation  velocity,  u/U^.,  since  the  agree¬ 
ment  in  the  z -component,  w/U(Jl,  is  generally  very  good  for  all 
Mach  numbers. 

In  Figure  8g  results  are  shown  for  a  very  thick  ellipse, 

£  =  0.8,  in  approximately  incompressible  flow,  Mco=0.1.  These 
results  show  that  the  first  order  method  is  good  regardless 


44 


of  thickness  for  very  Low  Mach  numbers.  (Again,  all  solu¬ 
tions  are  nearly  idea  Lira!  for-  approximately  incompressible 
flows,  being  exactly  identical  for  Mu,~0.  )  Finally,  in  F ip; 
Bh,  results  are  shown  for  a  thick  ellipse,  £=0.5,  at  a  free 
stream  Mach  number  of  0.5.  As  expected  the  comparison  is 
not  good,  since  the  method  was  developed  assuming  small  per¬ 
turbations  (slender  bodies)  for  the  outer  solution,  that 
being  the  dominant  solution  of  the  composite.  Note  however, 
the  inner  solution  remains  quite  acceptable. 

The  results  of  these  comparisons  demonstrate  that  the 
first  order  method  is  a  valid  approach  for  obtaining  accu¬ 
rate  and  uniformly  valid  solutions  about  slender  bodies  with 
round  edges  in  subsonic  compressible  flow.  The  first  order 
method  is  seen  to  be  easily  applicable  for  arbitrary  config¬ 
urations  using  numerical  techniques  (Chapter  V)  that  solve 
the  Prandtl-Glauert  equation,  with  the  mass  flux  surface 
boundary  condition.  The  second  order  application  for  the 
2-D  ellipse  is  presented  in  Appendix  D  while  surface  values 
of  the  term  ft  <j);i  as  obtained  by  various  approximations 
using  the  first  order  solutions  are  presented  in  Appendix 
E.  The  latter  comparisons  again  attest  to  the  validity  of 
the  first  order  method  in  obtaining  uniformly  valid  solu¬ 
tions  as  demonstrated  in  this  chapter. 
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2-D  Ellipse 


-  v/U^,  Thin  Airfoil  Theory  | 

I 

Jameson 
o  V/U0. 
a  tu(  /Uoj 

C  iwt  /U^ 

1st  Order  Multiplicative 
Composition  (A=l  )  ! 

-  V/lL  | 

-  (Ui/u*  J 

-  twl  /UM 


t 


46 
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Figure  8a.  Total  and  Perturbation  Surface  Velocity  Com¬ 
parisons,  First  Order  Analytic  Multiplicative  Composite 
(A=l),  Thin  Airfoil  Theory,  and  Jameson  Solutions,  2-D 
Ellipse,  0-50%  Chord,  €=0.1,  M^O.l 
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SURFACE  VELOCITY  (v/U^,  lu/U*  ,  iw/U, 
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Figure  8b.  (continued),  £=0.1,  1^  =  0.  3 
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SURFACE  VELOCITY  (V/U^,  lu/U,  ,  iw/U 


Figure  8c.  (continued),  £=0.1,  M  =0.5 


SURFACE  VELOCITY  (V/U^,  lu/U® ,  »w/U, 


Figure  Bd.  (continued),  £=0.1,  M  =0.6 
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Figure  8e.  (continued),  £  =  0.1,  Mro=0.7 
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SURFACE  VELOCITY  ( V/U^ ,  lu/U*  ,  <w/U£ 


Figure  8f ,  (continued),  £=0.1,  M  =0.8 


SURFACE  VELOCITY  ( V/t 


Figure  8h.  (continued),  M  =u.  S 


IV  Numerical  Application 
of  the 

F  i  r.'i  I.  Order  Method 


The  purpose  of  this  chapter  is  to  demonstrate  the  numer¬ 
ical  application  of  the  method  by  applying,  the  first  order 
method  numerically  to  both  2-D  and  3-D  bodies  in  subsonic 
compressible  flow.  These  results  are  obtained  at  zero  degrees 
angle  of  attack.  The  numerical  technique  used  was  a  Woodward 
code,  USSAERO  version  B  (Ref  2),  which  solves  the  Prandtl- 
Glauert  equation  using  distributed  surface  singularities. 

Both  the  velocity  surface  boundary  condition  (VBC)  and  the 
mass  flux  surface  boundary  condition  (MFBC)  were  applied  with 
the  latter  of  course  being  used  to  apply  the  first  order  meth¬ 
od.  First,  the  first  order  method  is  applied  to  a  NACA  0012 
airfoil  and  then  to  a  3-D  glide  bomb  configuration.  For  the 
airfoil  the  first  order  composite  results  are  compared  to  the 
Woodward  results  (VBC  and  MFBC),  Jameson  results,  thin  air¬ 
foil  theory  (Appendix  C),  and  experimental  results  (Ref  25, 
26).  In  addition,  the  Woodward  results  (VBC  and  MFBC)  will 
be  compared  showing  the  MFBC  results  to  be  the  better  and 
more  correct  of  the  two  solutions.  For  the  glide  bomb  the 
first  order  results  will  be  compared  only  to  the  Woodward 
results.  In  both  cases  the  surface  pressures  are  used  to 
evaluate  the  inviscid  drag  coefficients  with  the  result  be¬ 
ing  a  constant  improvement  in  value  for  all  Mach  numbers  us¬ 
ing  the  first  order  method  as  compared  to  the  Woodward  re¬ 
sults.  The  purpose  of  this  demonstration  is  to  show  that  the 
first  order  method  can  be  easily  applied  to  any  numerical 
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technique  that  solves  the  Frandtl-Glauert  equation,  ouch  as 
the  Woodward  code,  and  that  a  uniformly  valid  and  more  accu¬ 
rate  solution  results. 


Statement  of  Problem/Method  of  Solution 

The  first  order  governing  equation  for  2-D  problems  is 


RVL,  -o 


while  for  3-D  problems  it  is 


where  as  before 


f/3%  Outer  (A -o') 

fV'  -  j  A  i  iritermertiale.  (0<  A 2. 2  ) 

(_  l  j  Lnsuir  ( A  -2\ 


(IV-1) 


(IV-2) 


(lv-3) 


For  all  applications  to  follow,  A=1 ,  for  the  intermediate 
problem.  Using  the  Woodward  code,  standard  results  (outer 
or  Prandtl-Glauert  equation)  are  obtained  running  the  con¬ 
figuration  at  the  free  stream  Mach  number.  An  intermediate 
solution  for  the  same  configuration  is  obtained  for  an  inter¬ 


mediate  Mach  number,  Mt,  given  by 

Mc:/77A  (IV-4) 

for  a  given  value  of  A.  Finally,  the  inner  solution  is  ob¬ 
tained  for  the  same  configuration  in  incompressible  flow, 
M=0,  The  first  order  composite  solution  is  then  formed  us¬ 


ing  Eq  (11-39)  or 


and  the  surface  pressures  (for  all  solutions)  are  calculated 
using  the  isentropic  pressure  relationship, 

c,,  "DniO ■£)/'■' -/}  ov-5) 

Once  the  respective  solutions  have  been  obtained,  the  invis- 
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cid  drag  coefficients  will  be  evaluated  using  the  surface 


o 


pressuro  distributions.  For  the  Woodward  code,  the  accuracy 
of  these  calculations  of  course  depends  on  the  number  of  pan¬ 
els  or  data  points,  panel  density  and  distribution.  No  at¬ 
tempt  was  made  to  optimize  the  panel  density  and/or  distribu¬ 
tion  for  these  calculations  as  only  the  qualitative  compari¬ 
sons  are  of  interest. 


NACA  0012  Airfoil 

The  first  numerical  application  of  the  first  order  method 
was  to  a  NACA  0012  airfoil  (2-D)  whose  surface  is  given  by 
(Ref  24), 

-J1GX  -.351  x*  K2P4k3-.IC£x4}  (IV-6) 

where  Zt(x)  is  the  thickness  distribution.  The  configuration 
was  run  using  the  Woodward  code  (MFBC  and  VBC)  with  the  max¬ 
imum  number  of  chordwise  panels  (30)  at  free  stream  Mach  num¬ 
bers  of  0.4,  0.486,  0,6,  and  0.71,  all  at  zero  angle  of  attack. 


I 


Figure  9.  NACA  0012  Airfoil 
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Tim  MKMC  Woodwiii'ti  :n>lul  ion:;  I'orm  Llm  oulitr  or*  I'rnrul  L 1  -C  lauorL 
solution  used  in  the  f  omn  t.  i  on  of'  I. ho  composite}  solution.  The 
same  configuration  was  then  run  (MFBC)  at  intermediate  Mach 
numbers  of  0.198,  0.239,  0.292,  and  0.342  for  the  intermediate 
solution  (A=l),  and  then  finally  at  M=0  for  the  inner  solu¬ 
tion.  The  solutions  were  then  combined  to  form  uniformly 
valid  solutions  using  the  multiplicative  composition,  £q  (II— 
40),  in  terms  of  the  surface  velocity.  Results  were  compared 
to  the  Jameson  code,  analytic  thin  airfoil  theory  (Appendix 
C)  as  well  as  experimental  data. 

Results  of  the  numerical  application  of  the  first  order 
method  are  shown  in  Figs  10a  through  lOd.  Shown  are  the  sur¬ 
face  pressure  distributions  (Cp)  for  the  free  stream  Mach 
^  numbers  of  0.4,  0.486,  0.6,  and  0.71.  At  each  Mach  number 

the  composite  results  are  compared  to  the  VBC  and  KFBC  Wood¬ 
ward  results,  Jameson,  thin  airfoil  theory,  and  experimental 
free  flight  or  wind  tunnel  data  (Ref  25,26).  For  all  the 
Mach  numbers,  the  thin  airfoil  theory  results  are  clearly  in 
error  in  the  stagnation  region.  However,  these  results  are 
a  reasonable  approximation  to  the  adverse  pressure  gradient. 
The  Jameson  or  exact  numerical  solution,  as  compared  to  the 
Woodward  (VBC  and  MFBC )  and  composite  solutions,  is  seen  to 
have  higher  negative  peak  pressures  near  the  point  of  maxi¬ 
mum  surface  velocity  which  becomes  more  apparent;  as  Mach  num¬ 
ber  increases.  In  comparing  the  two  Woodward  solutions  it 
is  readily  seen  as  Mach  number  increases  that  the  MFBC  solu- 
tion  is  the  better  solution  with  respect  to  the  Jameson  re- 
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;;uLl:i.  noth  :»<>  Lu  t  i  uii:i  yieLd  a  pprox  i  mu  le  l.y  the  name  maximum 

r-' 

iiiirt'iii'i'  vulorl  l.y  tail.  Ilia  Ml-’IU!  nolulion  In  a  hoi. tor  match  t.o 
the  Jameson  solution  in  the  stagnation  region.  Thin  differ¬ 
ence  in  the  stagnation  region  will  show  up  dramatically  in 
the  inviscid  drag  calculations.  The  first  order  composite 
solutions,  being  formed  using  the  MPBC  Woodward  solutions, 
are  also  a  very  good  match  in  the  stagnation  region  to  the 
Jameson  solution,  especially  since  they  are  identical  at  the 
stagnation  point.  The  ability  to  correctly  model  the  flow 
in  the  stagnation  region,  as  previously  pointed  out,  is  im¬ 
portant  as  the  results  of  the  inviscid  drag  calculations  will 
demonstrate.  One  limitation  of  the  first  order  method  and 
Woodward  solutions  is  clearly  shown  in  Fig  lOd  in  the  inabil- 

^  ity  to  accurately  predict  the  maximum  surface  velocity  (Jame¬ 

son  solution)  as  the  flow  approaches  the  transonic  regime. 

As  previously  noted  these  first  order  methods  are  only  valid 
for  subsonic  or  supersonic  flows.  The  experimental  data  is 
matched  well  by  all  solutions  except  at  the  trailing  edge  in 
Figs  10b  and  10d  where  the  viscous  effects  become  apparent, 
and  only  the  Jameson  solution  is  a  good  match  to  this  data 
as  the  flow  approaches  transonic  Mach  numbers  (Fig  lOd). 

These  comparisons  verify  the  fact  that  inviscid  methods  are 
able  to  accurately  model  real  flow  behavior  when  used  pru¬ 
dently. 

The  resultant  inviscid  drag  coefficients  per  unit  span 
were  determined  using  the  pressure  distributions  shown  in 
Figs  10a  through  lOd,  the  theoretical  value  of  course  being 
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Figure  10a.  Comparison  of  Surface  Pressure  Distributions, 
First  Order  Multiplicative  Composite  (A=l),  Woodward  (VBC 
and  MFBC),  Jameson,  Thin  Airfoil  Theory  Solutions,  and 
Experimental  Data,  NACA  0012  Airfoil,  cx~ 0  deg,  MflP=0.4 
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Figure  10b.  (continued),  Mft)=0.48b 
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Figure  10c.  (continued),  M  =0,6 
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Table  1J 


Comparison  of  Drag  Coefficient.:;  per  unit  span 
NACA  U012  Airfoil 


r^ 


0.4 

0. 48o 

0.6 

0.71 

Jameson 

-.  000062 

-.000077 

-.  000109 

-.000162 

Woodward,  VBC 

.00513 

.00691 

.  01 1 36 

. OI656 

Woodward,  MFBC 
Cdw 

-.0022? 

-. 004o8 

-.01103 

-.02550 

1st  Order  (A=l  ) 
Composite , MFBC 
Cdc 

-.00005 

-.00158 

-.00555 

-.01408 

<iCd=  ICdw~Cdcl 

. 00222 

. 00310 

. 00548 

. 01142 

%  Reduction, 
ACd/ Cdw 

97.8 

06.2 

49.6 

44.7 

zero  for  2-D  inviscid  flow. 

These 

results  are 

shown  in  Table 

II  as  computed  for  the  Jameson,  Woodward  (VBC  and  MFBC ) ,  and 
composite  solutions,  where 


(IV-7) 


and 


D’-  Drag  per  unit  span  (lbf/ft) 

The  surprising  thing  to  note  in  Table  II  is  the  level  of  the 
inviscid  drag  coefficients  for  the  Woodward  results  in  light 
of  the  theoretical  result  being  zero.  With  these  levels  it 
is  easily  understood  why  difficulties  are  encountered  in  es¬ 
timating  inviscid  drag  levels  for  3-D  configurations  (drag 
due  to  lift).  Regardless,  it  is  seen  that  this  effort  can 
be  appreciably  aided  by  using  the  uniformly  valid  composite 
solution.  Although  the  differences  in  pressures  (Figs  10a 
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through  lOd)  between  the  Woodward  and  composite  solutions 
appear  negligible,  there  is  a  consistent  reduction  in  the 
drag  coefficients  for  all  Mach  numbers.  This  reduction,  in 
the  direction  of  the  exact  solution,  is  directly  due  to  the 
correction  made  in  the  stagnation  region,  thus  the  direct 
benefit  of  using  a  uniformly  valid  composite  solution  has 
been  demonstrated.  Being  able  to  better  predict  inviscid 
drag  levels  for  initial  performance  estimates  (including 
viscous  drag)  and  using  uniformly  valid  pressure  distribu¬ 
tions  only  enhances  the  design  process. 

Before  leaving  this  table  it  is  interesting  to  note  that 
when  the  VBG  is  applied  using  the  Woodward  code,  positive 
values  of  C^  result  while  when  using  the  MFBC  negative  values 
of  Cd  result.  Recall  that  in  Chapter  I  it  was  stated  that 
when  using  thin  airfoil  theory  results  to  calculate  the  in¬ 
viscid  drag  of  a  nonlifting  body  in  compressible  flow,  it  be¬ 
comes  necessary  to  add  a  leading  edge  drag  term  to  obtain  the 
correct  inviscid  result  of  zero.  In  other  words  the  thin 
airfoil  theory  solution  itself  results  in  a  thrusting  body 
which  is  consistent  with  the  results  obtained  using  the  Wood¬ 
ward  code  with  the  MFBC.  This  again  leads  one  to  believe 
that  that  the  MFBC  is  the  correct  equation  to  be  used  for 
these  approximations. 

Glide  Bomb 

The  final  application  of  the  first  order  method  is  to 
the  3-D  glide  bomb  configuration  (Ref  27)  shown  in  Fig  11. 
This  configuration  (originally  a  proposed  submunition)  is 


typical  of  the  current  glide  bomb  (and  cruise  missile)  con- 
I'igurit  ti  011:1  lluit  operate  in  the  subsume  speed  runge.  it  has 
a  hem  1 s pho r i ca l  nose  with  a  cylindrical  body  and  aft  boat- 
tail.  The  planar  w  i  ng ,  in  this  case  with  NAC.A  001  2  airfoil 
sections,  is  mounted  on  the  top  of  the  body.  The  cruciform 
tail  consists  of  four  flat  plates  with  round  leading  edges. 
The  configuration  was  run  at  a  zero  angle  of  attack  using  the 
Woodward  code  (MFBC)  at  free  stream  Mach  numbers  of  0.3»  0.4, 
0.5#  and  0.6.  The  corresponding  intermediate  solutions  (A=l  ) 
were  obtained  at  intermediate  Mach  numbers  of  0.149,  0.196, 
0.245,  and  0.292.  The  inner  solution  was  again  obtained  for 
the  incompressible  flow  case,  M=0.  Surface  pressure  data  are 
presented  in  Figures  12  and  13  for  the  Woodward  and  first 
order  multiplicative  composite  solutions  for  the  free  stream 
Mach  number  of  0.5.  These  data  are  representative  of  the 
data  at  the  other  Mach  numbers.  In  Fig  12  is  presented  the 
body  surface  pressures  and  Fig  13  presents  the  wing  upper 
and  lower  surface  pressures.  The  total  number  of  body  panels 
used  was  180  and  the  total  number  of  wing  and  tail  panels 
used  was  96.  No  attempt  was  made  to  increase  panel  density 
to  improve  the  results  as  would  occur  in  the  design  process. 

The  body  surface  pressures  shown  in  Fig  12  are  at  vari¬ 
ous  meridan  angles  about  the  body.  There  is  very  little  dif¬ 
ference  in  the  results  except  in  the  stagnation  regions  (both 
the  nose  of  the  body  and  near  the  wing  leading  edge),  as  be¬ 
fore  with  the  NACA  0012  airfoil  data.  The  same  is  true  of 
the  wing  surface  pressure  data  as  shown  in  Fig  13  for  three 


07 


I 


r* 


-1 . 0 


0.  0 


o°o 

o 


n 

o 

P 

o 


1.0 


cc 

°  -1.0 

w* 

ac 

to 

A 

M 

cc 

^  0.0 

6. 

o 

w 

M 

CD 

S  1*0, 

w 

o 

O  -1.0 


0 


0.  0 


©0 


c 

0 

n 


o 


i .  o 


,F 

< ) 


n _ D _ O 

0.5 


o 


Q 


x/1 


O —  0  =  15  deg 

1.0 


0°G 

0 

O 

G 

G 

0  Q _ 

_ Q 

_ G _ . 

G 

0.5 

x/1 

1. 

D 

0=45  deg 


Glide  Bomb 
0=0  deg,  M*=0.5 

O  1st  Order  Mult. 

Composite,  A=1 
□  Woodward,  MFBC 


Q  O  O 
0.5 


G 


x/1 


__0 — ,  0=75  deg 

1. 0 


Figure  12.  Comparison  of  Body  Surface  Pressure  Distri¬ 
butions,  First  Order  Multiplicative  Composite  (A=l  )  and 
Woodward  (MFBC)  Solutions,  Glide  Bomb,  o  =  0  deg,  ^=0.5 
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COEFFICIENT  OF  PRESSURE, 


Figure  12.  (continued) 
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Figure  13.  Comparison  of  Wing  Upper  and  Lower  Surface 
Pressure  Distributions,  First  Order  Multiplicative 
Composite  (A  =  l )  and  Woodward  (MFBC)  Solutions,  Glide 
Bomb,  cx= 0  deg,  1^=0.  5 
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spanwise  stations  on  the  right  wing.  For  the  body  data  the 
tiaino  t.rmuh;  are  observed  nn  before  for  the  flow  about  flu; 
nose  or  stagnation  region.  The  composite  results  show  a 
slightly  higher  velocity  in  this  region  and  reduced  negative 
peak  pressures  as  compared  to  the  Woodward  results.  For  the 
wing  data,  Fig  13,  "the  velocities  in  the  stagnation  region 
for  the  composite  solution  are  again  slightly  higher  than  the 
Woodward  results,  however,  the  inboard  station  shows  a  slight¬ 
ly  greater  negative  peak  pressure  while  the  outboard  station 
is  slightly  less  than  the  Woodward  results.  One  can  clearly 
see  in  Fig  12  the  ever  increasing  presence  of  the  wing  as  one 
moves  around  the  body,  while  in  Fig  13  the  outboard  loading 
is  somewhat  less  than  the  inboard  due  to  the  presence  of  the 
wing  tip. 

The  resultant  drag  coefficients  were  calculated  using 
the  surface  pressure  data  where 

fiDls  <IV'8> 

and 


D  -  Pressure  or  Form  Drag  (lbf) 
with  the  results  being  presented  in  Table  III.  Although  the 
conf iguration  was  at  a  zero  angle  of  attack  a  negligible  lift 
was  produced  due  to  the  high  mounted  wing  and  resulting  wing 
body  interference.  As  shown  in  Table  III  there  is  again  a 
consistent  reduction  in  the  value  of  the  drag  coefficients 
even  though  the  difference  in  pressure  distributions  was  a- 
gain  slight.  This  is  directly  the  result  of  using  a  solution 
(composite)  which  more  accurately  models  the  flow  in  stagna- 
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Comparison  of  Drag  Coefficients 
Glide  Bomb  Configuration 


0.3 

0.4 

2-1 

0.  6 

Woodward,  MPBC 
°Dw 

-.06822 

-.07026 

-.07449 

-.08364 

1st  Order  (A  =  l ) 
Composite, MFBC 

cDc 

-.06736 

-.06869 

-.07252 

-.07917 

lACpI.  Drag 

Counts 

8.6 

15.7 

19.7 

44.7 

tion  regions.  Although  the  difference  in  pressure  distribu¬ 
tions  is  slight,  these  differences  are  important  for  drag 
calculations.  Note  again  the  relatively  high  values  for  the 
inviscid  drag  coefficients  in  light  of  the  theoretical  value 
being  zero  for  nonlifting  bodies.  The  results  show  the  first 
order  method  to  be  an  easily  applicable  and  accurate  method 
for  obtaining  uniformly  valid  solutions  about  2-D  and  3-D  con¬ 
figurations  with  round  leading/trailing  edges  in  subsonic  po¬ 
tential  flow.  The  proposed  method  has  also  been  demonstrated 
to  be  easily  applicable  to  a  current  numerical  technique  solv¬ 
ing  the  Prandtl-Glauert  equation.  Application  of  the  proposed 
method  has  resulted  in  consistently  improved  inviscid  drag 
estimates.  Although  the  results  presented  have  been  for  non¬ 
lifting  bodies  only,  the  first  order  method  can  be  easily  ex¬ 
tended  to  lifting  problems  (camber  and  or  angle  of  attack), 
however  this  remains  to  be  verified. 


V  Conclusion:: 


The  purpose  of  this  dissertation  was  to  develop  a  method, 
based  on  the  principles  of  singular  perturbation  theory,  for 
obtaining  uniformly  valid  solutions  about  slender  bodies  with 
round  leading/trailing  edges  in  subsonic  potential  flow.  The 
first  order  method  an  developed  using  the  multiplicative  com¬ 
position  has  been  shown  to  satisfy  this  goal.  Specifically: 

1.  The  first  order  method  corrects  small  perturbation  theory 
(outer  or  Prandtl-Glauert  solution)  with  a  stagnation  region 
theory  (inner  solution)  yielding  uniformly  valid  solutions. 

The  method  was  developed  by  formulating  approximate  equations 
to  the  exact  potential  equation,  and  mass  flux  surface  and 
infinity  boundary  conditions.  The  mass  flux  surface  boundary 
condition  was  shown  to  be  the  boundary  condition  to  use  for 
these  approximations. 

2.  Analytic  solutions  were  developed  for  application  of  the 
first  order  method  to  the  2-D  ellipse  with  the  results  being 
compared  with  thin  airfoil  theory  and  a  numerical  solution 

of  the  exact  potential  equation  as  developed  by  Jameson.  Com¬ 
parisons  with  the  Jameson  results  were  very  favorable  within 
the  limitations  of  the  theory,  that  is  for  slender  bodies  in 
subsonic  compressible  flow. 

3.  The  first  order  method  was  then  applied  numerically,  using 
a  panel  method  as  developed  by  Woodward,  to  a  NACA  0012  air¬ 
foil  with  the  results  being  compared  again  to  Jameson  results. 
Woodward  results,  thin  airfoil  theory,  and  experimental  data. 
The  comparisons  again  proved  to  be  very  favorable  within  the 


limitations  of  the  theory.  Next,  the  first  order  method  was 
applied  numerically  to  a  '3-D  glide  bomb  configuration  using 
the  Woodward  code  with  the  results  being  compared  to  the  Wood¬ 
ward  or  outer  solution  results.  The  numerical  applications 
of  the  first  order  method  demonstrated  the  ease  of  applica¬ 
bility  and  the  versatility  of  this  proposed  method  for  any 
numerical  technique  solving  the  Prandtl-Glauert  equation  with 
the  appropriate  boundary  conditions. 

4.  The  latter  two  numerical  applications  of  the  first  order 
method  resulted  in  consistent  improvements  in  the  prediction 
of  inviscid  drag  levels  as  compared  to  the  Woodward  results. 
This  increase  in  accuracy  leads  directly  to  increased  confi¬ 
dence  in  inviscid  aerodynamic  predictions  and  their  later  use 
when  predicting  viscous  effects.  This  increase  in  accuracy 
has  been  shown  to  be  directly  due  to  the  correction  made  to 
the  outer  solution  in  the  stagnation  region  by  using  the  com¬ 
posite  solution. 

5.  Finally,  and  most  importantly,  the  method  as  proposed 
obviates  the  need  for  analytic  solutions  or  numerical  patch¬ 
ing  of  solutions  in  forming  a  uniformly  valid  composite  solu¬ 
tion.  Three  solutions  are  required  for  the  composite  solu¬ 
tion.  all  being  solutions  to  the  Prandtl-Glauert  equation  for 
the  same  surface  at  three  different  Mach  numbers  (M^,  Mt,  and 
M-0).  The  composite  solution  once  obtained  offers  increased 
accuracy  of  simulation  as  demonstrated  in  this  effort. 

Before  closing  one  item  should  be  remembered  in  light  of 
the  results.  This  concerns  the  perturbation  parameter,  £, 
which  was  assumed  to  be  small  throughout  this  effort.  Small 


is  a  relative  term,  however,  usual  definitions  (Ref  7,19)  im¬ 
ply  that  the  perturbation  parameter  is  on  the  order  of  ten 
percent  (or  less)  of  the  parameter  being  perturbed.  For  prac¬ 
tical  application,  however,  an  upper  limit  of  fifteen  percent 
is  often  used.  Applications  beyond  these  limits  are  at  the 
expense  of  reduced  accuracy.  All  applications  of  the  proposed 
method  in  this  effort  have  been  for  thickness  ratios  of  ten 
percent  or  greater  and  thus  the  method  has  been  tested  for  the 
extreme  cases  in  terms  of  the  perturbation  parameter  and  accu¬ 
racy.  For  application  of  the  method  to  thinner  bodies  the 
results  should  be  even  better. 

To  summarize  then,  the  first  order  method  is  a  highly 
accurate  and  easily  applicable  technique  for  obtaining  uni¬ 
formly  valid  results  about  slender  bodies  with  round  leading/ 
trailing  edges  in  compressible  subsonic  potential  flow.  The 
method  as  developed  was  applied  to  nonlifting  bodies  only, 
however,  the  method  can  and  should  be  developed  for  lifting 
bodies  (camber  and  or  angle  of  attack). 


n 
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Appendix  A:  First  Order  Analytic 
Solution  for  a  n  Jill,  ipt  1  c  Cylinder 
in  Subsonic  Flow 


As  presented  in  Chapter  III,  Eqs  ( III— 2 )  through  ( 1 1 1 - 5 ) , 
the  outer,  inner,  and  intermediate  first  order  problems  for 
the  elliptic  cylinder  (Fig  2),  whose  surface  is  given  by 

x '  >  m- 1  =  o  ( a-i ) 

in  the  subsonic  compressible  isentropic  flow  of  a  perfect  gas 
(Y=1.4)  are  all  of  the  form 


+  0tLi  -  O  (A-2) 

on  F(x,z)  =  0  (A-3) 

-■>  O  at  infinity  (A-4) 

where 

(fil,  outer  (K  - c)  .  i  --  4 


r  \A  ,  ulte.rwecLute  (A-5) 

[/  ,  i AAe,-  ( h  o)  ,  2' -6’1 

To  transform  the  problem  to  the  incompressible  flow  case  a- 
bout  an  equivalent  elliptic  cylinder  let  k-Yg,,  ip where 
B*0),  and  <$>■?(/-  %Bl)x  resulting  in 

-  O  (A-6) 

$,i  on  F(x,  z  )=0  (A-7) 

K  Yj  at  infinity  (A-8) 

The  problem  is  now  mapped  to  a  known  solution,  the  circular 
cylinder  in  incompressible  flow,  using  the  Joukowski  trans¬ 
form  and  its  inverse.  First,  however,  it  is  necessary  to 
determine  the  Joukowski  transform  constants  when  mapping  from 
the  circle  to  the  ellipse  where 
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The  known  solution,  represented  by  the  complex  potential,  in 
the  plane  is 

J /?)'-  $//}»<  i-’U) 


(A-9) 

The  circle  radius,  A »  and  the  constant  C  in  the 

Joukowski 

transform  are  determined  to  be 

fi = ie ( 1  * 

(A-10) 

c  uyi  -s'td 

Next,  the  inverse  transform  is  determined  to  be 

( A-ll ) 

(A-12) 

It  becomes  necessary  at  this  point  to  introduce 

elliptic  coor 

dinates  so  that  the  inverse  transform  and  resulting  complex 

potential  may  be  separated  into  distinct  real  and  imaginary 

parts  (this  cannot  be  done  in  Eq  (A-12)).  Let 

V  =  kch  /  ,  V  -  kdlK' 

(A-13) 

where 

(A-14) 

resulting  in 

/  :  k  cfyu  ‘  C'oV  vJhexe  ch  ?  cc-^lx 

(A-l 5 ) 

H  -  k  V  inhere  sk  f  smA 

(A-16) 
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O  v<tjn  ,  >/  *  jl(  (A-17  ) 


It  is  interesting  to  note  that  in  elliptic  coordinates  the 
surface  boundary  condition  is  simply 


(A-18) 

where  the  total 

velocity  is  given  by 

Since 

v  1  1 

(A-19) 

and  the  surface 

--  /  --  (\~j-  )  * 

of  the  ellipse  is  given  by 

(A-20) 

(6;V<-  (%■)'  --  i 

then 

<-yk.  --Vnk'  Aa  = 6 A 

(A-21  ) 

or 

specifying  the 

n  .  /'T&f 

/'<o  -  ..fr?  1/  ,  _  ^ 

surface  F(x,z)=0.  From 

(A-22 ) 

-■  ; 

( A— 2  3 ) 

we  find 

t  -  ;v  fl  ■  (frV 

-  f4C  ’ 

(A-24) 

From 

y/./Z-i  P  C./t  1  [,  k  (*  *  1  *■  ^  l 

( A — 2  5  ) 
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we  find 


/( -  d)  1  k) 

V  -  CoS' '(  ) 


where 


(A-20 ) 
(A-27) 


g  .  1/  2  J  >  ^  v  1  I, )  ' 

£  ■  f/1  •  (s  L)‘ 

Finally,  the  inverse  transform  is  given  by 

J(l-)  r  4  er 

The  potential  is  then  given  by 

&/>.*')-  I  le"  *■  (  *■ 

Note,  the  arbitrary  constant  in  Eq  (A-31 )  results  since  this 
is  a  Neumann  problem.  Returning  to  the  (x,z)  coordinates, 
the  first  order  analytic  solution  for  an  elliptic  cylinder 
is  given  by 

0  (*,?.)  -  )f;  [<?+{¥^r]h*Rl]  -  (7-  y)|.  *■  (A-32) 


(A-28) 

(A-29) 

(A-30) 

(A-31) 


where 


i-;. p  ■(/,’( 


( A  -  3  3 ) 

(A-34) 
(A-35 ) 
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Appendix  B:  Thin  Airl’oil  Theory 
Solution  for  an  Elliptic  Cylinder 
in  Subsonic  Flow 


The  equations  for  the  flow  over  an  elliptic  cylinder 
(Fig  2)  in  incompressible  flow  given  by  thin  airfoil  theory 


'/L  f'  'V y f  --o  (B-i) 

-  04  (B-2) 

Ucf]'-,  C)  <tt  (B-3) 

Assuming  the  solution  can  be  given  by  the  following  expres¬ 


sion, 


cpu,i\  --  ~  [q( $-,c)  l  fiHYl  I'cli  (B-4) 
which  represents  a  source  distribution  on  the  x-axis,  /x/^1, 
where  Q(^',0)  is  the  source  strength.  The  perturbation  com¬ 


ponents  of  velocity  are  given  by t 


,r  K  ±  ['  0(y.c)(x  -?_)  /, 
((u  '  ,377 J,  ITfVTk*  •’ 

,a  _  W.  ..If  _£jQl(l<zI  c-h 

"  Uu,~  -hi  j.,  la-jV  ^  J  1 


(B-5) 

(B-b) 


To  first  order  from  conservation  of  mass  the  source  strength 


is  given  by: 


Q  (*.  c)  -  ^  ^ 


(B-7) 


The  surface  perturbation  velocity  is  thus  given  by, 


,  X  f"  rBi'-c/. 
1  '  J;  (  *  -  1 


ij 


r  f  / 


( B-8 ) 


with  the  surface  velocity  being: 

il  r  !  (B_9) 

Applying  the  Prandtl-Glauert  compressibility  correction  gives 
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the  following  subsonic  surface  velocity j 

rl ;  i ac'  ‘b-10> 

It  is  quickly  noted  that  this  value  is  constant  for  a  given 
free  stream  Mach  number, 
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Appendix  Ci  Thin  Airfoil  Theory 
Solution  for  n  NASA  001?  Airfoil 
in  Subsonic  Flow 


Following  the  procedure  of  Appendix  B,  the  surface  per¬ 
turbation  velocity  in  incompressible  flow  is  g iven  byi 

«  ^  -  i  -  i:-  (,  ti  t  •  &  /,!■(*  •»  </j  (c-'  » 

The  surface  of  the  airfoil  (Pig  9)  is  described  by, 

It-  -  j(,2-X2 /xTT -JAG  /p,7x  -.dJbx*  f  6U  y3 -.0&8x*J  (C-2) 

where  / x( -  1.  Applying  the  surface  boundary  condition,  Eq 


(B-7),  results  ins 

TLr  nUi  '^ftfrru-s))  - 

(<-y>  i  <  -}]  y1* 

*  *{<■•■*&•  >4-M'  * 

(  x  'f-t  -M-.  i  6.3 U  -  .032 <")  t-  (■£ )(  CC6 )  J 


(C-3 ) 


The  resulting  surface  velocity  is, 

£■■!>■  1  (c-4) 

with  the  result  for  compressible  flow  being  given  by s 

■V  11  Vwjr-Hi  (c-5) 

This  expression  for  the  surface  velocity  results  in  infinite 
or  undefined  velocities  at  the  leading  and  trailing  edges  of 
the  given  airfoil. 
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Appendix  Di  Second  Under  Analytic 
Solution  for  a n  Elliptic  Cylinder 
In  Subsonic  Flow 


r' 


For  application  of  the  second  order  method,  consider  the 
nonlifting  elliptic  cylinder  (Pig  2)  in  the  subsonic  compress¬ 
ible  isentropic  flow  of  a  perfect  gas  (Y=1.4).  The  surface 
of  the  ellipse  is  given  by  Eq  (III— 1  ).  The  second  order  out¬ 
er  ( A=0 ) ,  intermediate  ((KA<2),  and  inner  problems  are  of  the 
form 


where 


"cvv/:.. 

( D-l  ) 

(W',u)x  f  'fji  =o 

on  P(x,z)=0 

(D-2) 

/  tfAl  ->  0 

at  infinity 

(D-3) 

f  AC  ,  (-'“ter  U\  -c) 

/)/  -  j  JS  f  ,  l  r]  lifi’r'KC  lu?  Le 

r( 

V 

V 
O 

(D-4) 

O  ,  t. ')  ltd'  (  A  -■  J ) 


X  -  ( ■/hi)  -  ( Y- 1)  Sl  (D-5) 

It  can  be  shown  that  the  solution  to  this  problem,  if  it 
exists,  is  unique  to  within  an  unknown  constant  (e.g.,  let 
Y,  and  1/t  be  two  solutions  theny,  where  ^  - y,%  ,  is  equal 

to  a  constant).  From  the  Divergence  Theorem  for  infinite 
domains  (Ref  28)  one  can  also  show  the  following  compatibility 
condition  must  be  satisfied  for  a  solution  to  exist  (the  do¬ 
main  J2  representing  the  region  from  the  surface  of  the  ellipse 
to  some  circular  boundary  at  infinity)* 

j„  ( <■  II  =  O  (D-6) 

It  is  quickly  noted  that  this  condition  is  satisfied  for  a 
free  stream  Mach  number  of  zero,  or  incompressible  flow. 
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This  is  indeed  true  for  the  inner  approximation  where  the 
solution  can  easily  be  shown  to  be 


-  CoA*>Lint  (D-7) 

and  thus 

(O'8) 

Equation  (D-6)  represents  a  necessary  and  sufficient  condi¬ 
tion  for  the  existance  of  second  order  solutions.  For  com¬ 
pressible  flow  Eqn  (D-6)  can  be  shown  to  be  satisfied  for 
the  elliptic  cylinder  and  thus  we  are  assured  of  the  exist¬ 
ence  of  an  analytic  solution  for  the  second  order  outer  and 
intermediate  problems  and  a  second  order  composite  solution. 
This  in  turn  reinforces  the  credibility  of  the  first  order 
method  as  will  also  be  demonstrated  in  Appendix  E.  Given 
the  demonstrated  accuracy  of  the  first  order  method  however, 
plus  the  ease  of  obtaining  first  order  solutions  as  compared 
to  a  second  order  effort,  the  first  order  method  should  be 
sufficient  for  most  applications. 


Appendix  E:  Surface  Values  of  Bz al <+(£,*■ 
for  the  2-p  Ellipse  f 
in  Subsonic  Flow 

In  this  Appendix  some  interesting  comparisons  of  the  sur¬ 
face  values  of  the  term  /S2^^3f  as  evaluated  exactly  (numeri¬ 
cally)  and  by  various  approximations  are  shown  in  Figs  Ela-f 
(the  legend  for  these  Figs  is  in  Table  IV).  The  exact  value 
of  /6?<PKX+<paa  ,  obtained  numerically  using  results  from  the 
Jameson  code,  represents  the  exact  perturbation  potential 
equation,  Eq  (II-9).  The  approximations  are  developed  using 
the  first  order  results  with  being  formed  using  the 

additive  composition,  or, 

(tfU.i (fftx,?)  (E-i) 

Note  that  this  assumed  composite  solution  is  not  the  same  as 
the  composite  solutions  used  previously  for  the  velocity  (Eqs 
(11-39),  (11-40) ),  although  the  form  is  the  same.  Analytic 
results,  representing  various  approximations,  are  shown  for 
the  first  order  outer  or  Prandtl-Glauert  equation, 

(e-2  ) 

the  first  order  additive  composite  (A  =1 ) 

<E-3> 

and  the  second  order  additive  composite  (A=l )  using  the  first 
order  solutions 

PtL,  >C--  U/*<) -a-isp Jc.< 

05-4) 

As  seen  in  Figs  Ela-f  the  Prandtl-Glauert  equation  is  of 
course  equal  to  zero  for  all  Mach  numbers,  but  it  (as  well  as 


8) 


Table  IV 


Legend  for  Figures  Ela-f 
2-D  Ellipse 

o  Jameson , {$ t-  t 

-  1st  Order  Additive  Composite  (A=l ) ,  y3z4^  f  (fitc 

-  1st  Order  Additive  Composite  (A=l),  Approximation 

to  2nd  Order  Equation, 

-  1st  Order  Outer  Equation  =  O 

the  other  approximations)  is  an  excellent  approximation  com¬ 
pared  to  the  exact  solution  to  the  flow  over  the  majority  of 
the  body  where  the  assumption  of  small  perturbations  is  valid. 
The  exception  to  this  occurs  in  or  near  the  stagnation  region, 
and,  as  represented  by  the  exact  solution,  the  flowfield  is 
seen  to  be  distinctly  divided  into  two  regions.  The  first  and 
second  order  approximations  more  closely  approximate  the  exact 
solution  (especially  in  the  stagnation  region),  with  the  sec¬ 
ond  order  approximation  being  almost  identical  to  the  exact 
solution  for  all  Mach  numbers  except  at  1^=0. 8  (transonic  flow). 
The  change  in  character  of  the  exact  solution  as  one  approaches 
the  transonic  flow  regime  ( M^=  0.8)  clearly  shows  another  ap¬ 
proach  is  required  to  obtain  more  accurate  results  in  this  flow 
regime.  To  conclude  however,  these  comparisons  once  again 
demonstrate  that  the  method  as  developed  yields  solutions  that 
are  uniformly  valid,  however  it  is  applied. 


Figure  Ela.  Comparison  of  Surface  Values,  ,  Jameson,  First  Order  Outer 

(Thin  Airfoil  Theory),  First  Order  Additive  (A=l)  Composite,  and  First  Order 
Additive  (A=l )  Composite  Approximation  of  the  Second  Order  Equation,  2-D  Ellipse. 
0-50%  Chord,  £=0.1,  M  =0.1 
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